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ABSTRACT 


The Ocean Thermal Energy Conversion (OTEC) uses the low 
thermal energy potential available from ocean temperature 
gradients. A method is presented to analyze such systems and, 
for this purpose, a comprehensive simulation is developed. 

The simulation includes parasitic power requirements, losses 
due to interconnecting lines, and heat exchanger pressure 
drops. Cost functions are included and numerical optimization 
is employed to obtain optimal designs based upon minimum cost. 
The analysis is converted to a computer code and coupled to 
the COPES/CONMIN optimization code to facilitate a fully- 
automated design where the computer makes the design decisions 
and performance trade-off studies. The final product is an 
optimum power system module design for the designated net 
electrical output required and the specified system and design 
constraints. 

Preliminary results are presented for a range of system 
power levels. Optimum designs are obtained and compared for 
Systems in which either titanium or aluminum tubes are used 


in the heat exchangers. 
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I. INTRODUCTION 


A. BACKGROUND 

Qcean Thermal Energy Conversion (OTEC) is a concept 
using the low thermal energy potential available from the 
Ocean temperature gradient that exists between warm surface 
Ocean water and cold water in deep ocean regions. 

; The idea of converting the stored ocean energy to useful 
power originated with French physicist Jacques d'Arsonval in 
1881 [Ref. 1]. It was nearly a half-century later that the 
technical feasibility of ocean thermal energy conversion 
could be demonstrated. In 1926, George Claude used an open 
cycle power system to extract heat from surface water for 
indirect conversion of the thermal energy of a working fluid. 
Operating at a low pressure the working fluid was used to 
drive a turbine providing electrical power generation. 

Though Claude's limited power system produced only 
22 kilowatts of electricity while requiring approximately 
80 kilowatts of power to drive its equipment, it stirred 
the scientific and research community to consider the 
attractiveness of ocean thermal energy conversion [Ref. 2]. 
Claude called for immediate action on his ocean thermal 
power system, because of the Federal O11 Conservation Board's 
dire predictions that the United States had only six years 
of oil production remaining. Obviously the dire predictions 


ascribed to by the Federal 0il Conservation Board did not 


Mea 





come true, but the oil crisis of that period heightened 
scientific interest in extracting energy from the ocean. 

Now, 55 years later, the United States is faced with an 
energy crisis because of increasing industrial and social 
dependence on foreign petroleum. Dwindling supplies and 
erratic price hikes have rekindled interest in ocean thermal 
energy conversion, since it utilizes an inexhaustible supply 
of fuel. 

Currently, the United States Department of Energy is 
attempting to develop the necessary technology and demon- 
Strate the feasibility of large-scale OTEC power systems. 
However, there are major engineering development problems 
which must be solved before OTEC can be standardized and 
become a viable source of electrical power generation. 

The single controlling factor which creates troublesome 
technical encounters is low thermal power system efficiency 
(one to four percent depending upon parasitic power require- 
ments). Because the heat energy used by OTEC must be extracted 
from a small ocean temperature difference, extremely large 
volumes of surface water must pass through a proportionately 
Sized evaporator to provide sufficient indirect heat energy 
to convert the working fluid into vapor to drive a turbine- 
generator for electrical power generation. Concurrently, to 
convert the turbine exhaust to a saturated liquid, completing 
the closed cycle, a condenser having compatible heat absorp- 


tion capacity must be employed. 
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Economic handling of the volume of fluids required for the 
heat absorption, expansion, and heat rejection phases of the 
cycle requires close scrutiny of evaporator, turbine, condenser, 
and pump design to Minimize the parasitic losses with respect 
to the generated electrical output. Because of the low 
thermal efficiency, relative to nuclear or fossil fuel-fired 
power plants, the margins for design and operating error in 
QTEC plants will be narrow. 

With the advent of high-speed computers, numerical 
methods for solving these complex engineering problems with 
multiple design variables and constraints are now possible, 
The case for utilizing an optimizing scheme for not just 
one syStem component, but rather the complete power generation 
Gyele, can easily be made. In effect, it would serve as a 
Systems analysis tool, to optimize component design and cost, 
relative to a specific electrical output or to enable compari- 


sOn and evaluation of competing OTEC designs. 


Ba OssiECTIVES 

The objectives of this work are to develop a computer 
code for the Ocean Thermal Energy Conversion (OTEC) power 
System and to couple the analysis to a numerical optimization 
code to provide an optimum system design capability, considering 
both performance and economics. 

This would create an optimum modular design relative to 
a specified objective function for a desired net electrical 


output, such as a 25 MW (net) power system. Such a design 


IES 





would permit construction of higher capacity power systems 

using the optimized modules as substations of the total power 
plant. Cost savings, improved plant performance, redundancy, 
and reliability could be the immediate beneficiaries of such 


ae venture. 


See OVERVIEW OF THE OTEC POWER SYSTEM ANALYSIS 

To analyze the closed-cycle OTEC power system, the 
fundamental relationships of heat transfer, fluid mechanics 
and thermodynamics are used to simulate a variety of system 
component designs, which form the basis of the power system 
algorithm. The scope of this analysis will be limited to 
the OTEC power system and sea water systems only. Mooring 
Systems, power delivery, hull, and cold pipe design will not 
be addressed. 

The performance analysis will be divided into four 
sequential sections as shown in Figure 1, and discussed 
in detail in subsequent chapters of this thesis. 

Input parameters (design constants) for the power cycle 
analysis will include: 

Required net electrical output. 


salt water inlet temperature to the evaporator 
and condenser. 


Length of hot and cold salt water pipes. 
Heat exchanger tubing material (aluminum or titanium). 
Heat exchanger tube orientation and profile. 


Pump mechanical and motor efficiencies. 
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Turbine mechanical efficiency. 
Conoeteeneiloedantcal and electrical efficiency. 
bmormoulineg contre! factor. 

Piping absolute roughness. 


Projected annual inflation rate for aluminum heat 
exchanger retubing. 


is 





EVAPORATOR 


PERFORMANCE 


PARASITIC 


LOSSES 


TURBINE AND 


ELECTRICAL POWER 


REQUIREMENTS 


CONDENSER 


PERFORMANCE 





Figure 1. Power System Sequential Analysis 
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Peweonwenec ele DESCRIPTIONS 


A. INTRODUCTION 

This chapter will provide a brief description of the OTEC 
power system. First, looking at the ideal Rankine cycle, 
the fundamental thermodynamic concepts will be enumerated. 
Then the deviations from the ideal cycle will be presented, 
creating the configuration assumed for the present cycle 
analysis which will be amplified in detail by follow-on 


chapters. 


meee WeAt OTEC RANKINE CYCLE 

The closed-cycle OTEC concept is based upon a Rankine 
power cycle that is driven by the low thermal energy potential 
available from the ocean temperature gradient that exists 
between warm surface water and cold deep water in ocean 
regions. The power cycle consists of a working fluid circu- 
lation pump, evaporator (heat absorption), turbine (expansion), 
and condenser (heat rejection), as shown in Figure 2. The 
majority of current OTEC designs are based upon ammonia as 
the working fluid -- a design decision that is adopted for 
this analysis. 

Figure 2 also illustrates an ideal OTEC Rankine cycle, 
plotted on temperature-entropy coordinates. In the ideal 
cycle, the low pressure working fluid (state point 1) is 
i1sentropically pumped to the evaporator operating pressure 


(state point 2). The working fluid (ammonia) is then 


1 





NH; CIRC PUMP 


| 3 
EVAPORATOR 
: 
2 2 
4 
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Deere 2. lIdealized OTEC Rankine Cycle 
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SOmiemecad tO 4a Sdturated Vapor in the evaporator by indirect 
heat energy exchange from warm surface ocean water (state 
point 3). Mechanical power is generated by isentropic expan- 
Sion of the saturated ammonia vapor through the turbine 
(state point 4). 

After exiting the turbine, the wet, low-pressure vapor is 
converted to a saturated liquid in the condenser by indirect 
heat absorption from cold ocean water (state point 1), 


Metuiming the cycle back to the working fluid circulation 


pump. 


See aclUAL OTEC RANKINE CYCLE 
In actuality there are numerous deviations from the ideal 
cycle which must be considered in this analysis. These are: 
(1) Turbine, generator and pump efficiencies. 


(2) Pressure drops in evaporator and condenser (tube- 
Side and shellside). 


(3) Pressure drop across moisture Separator. 
(4) Elevation change and frictional losses in 
Piping: (a) re-flux pump piping, (b) piping 
from circulation pump to evaporator. 
(5) Evaperager outlet quality (85 to 95%). 
(6) Moisture separator outlet quality (99 to 99.5%). 
The deviations from the ideal Rankine cycle described 
above are depicted in the flow diagram and temperature-entropy 
plot of Figure 3. In the actual OTEC Rankine cycle, the low 
pressure working fluid (state point 1) is pumped up to the 


evaporator operating pressure by the ammonia circulation pump 


With an adiabatic efficiency (state point 2). The working 


1S, 





MOISTURE 







_ : SEPARATOR 
HOT New 
WATER 
RE-FLUX PUMP 
2 


NH , CIRC PUMP 


Bueuge 5. Nerual OFEC Rankine Cycle 
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fluid (ammonia) aren converted to a wet vapor with an 
evaporator outlet quality (85-95%) acting under a shellside 
pressure drop (state point 3). Evaporator outlet vapor then 
passes through a moisture Separator to improve vapor quality 
(99-99.5%) creating a pressure drop (state point 4). 
Mechanical power is generated by the expansion of the moisture 
Separator outlet vapor through the turbine with an adiabatic 
efficiency (state point 5). After exiting the turbine, the 
Wet, low pressure vapor is converted to a saturated liquid 
in the condenser acting under a shellside pressure drop 
(state point 1), returning the cycle to the working fluid 
circulation pump. 

This figure forms the thermodynamic basis for the OTEC 


power system analysis which follows, 


a: 





III. EVAPORATOR AND MOISTURE SEPARATOR 


A. INTRODUCTION 
several heat exchanger concepts have been proposed for 
Closed-cycle OTE€ systems. Among these designs are: 
Conventional shell and tube heat exchanger. 
Plate type heat exchanger. 
Within these basic concepts, variations in design have 
been proposed, including: 
Orienu@eaonmet tibes (horizontal or vertical ). 


Heat exchanger tube material (i.e., titanium, 
aluminum). 


. Method of tube enhancement (i.e., fluted, porous 
coatings). 


Location of tube enhancement (i.e., internal and/ 
or external). 


Location of the vapor separator (i.e., internal 
OY external). 


Location of the heat exchangers relative to the 
Sea surface. 


Method of biofouling control. 
The analysis. to be presented for the evaporative heat 
exchanger will be based on the following design characteristics: 
Single-pass shell and tube heat exchanger. 


Internal vapor separator with a gravity drain to 
evaporator inlet. 


Heorazontal Orientation of tubes with an equi- 
lateral triangle or square tube profile. 


Smooth plain-tube configuration (no enhancements). 


fay 





Tube material (titanium or aluminum based on a 
30-year life-cycle criterion). 


Biofouling control based upon an achievable 
POrmTrS factor - 


Heat exchanger centerline located on sea surface. 
As an overview of the evaporator-moisture separator 
analysis, the following major steps in the algorithm are 
listed in order of their execution (numbers in parentheses 
refer to equations developed in the subsequent analysis): 
eoeutileaeron Of System constants (see I.C.). 
Initialization of design variables (D.V.). 
Tube length. 
SW velocity through hot pipe. 
Inner diameter of hot pipe. 
Tube outer diameter. 
SW velocity through evaporator tubes. 
Inner diameter of NH3 piping. 
Inner diameter of NH3 re-flux piping. 
Tube profile pitch ratio. | 
plkpewater Mass flow rate (1) . 
Total number of tubes (2). 
Total heat transfer surface area (3). 


Assume an initial salt water bulk temperature (6), and 
ammonia heat transfer coefficient (9). 


Overall heat transfer coefficient (4). 
Number of transfer units (11). 
Heat exchanger effectiveness (13). 


Salt water outlet temperature (15). 


eS) 





Revised bulk temperature (16); iterate with (6). 
Amount of heat absorption (17). 

Log mean temperature difference (18). 

Film temperature (19). 


Initial ammonia mass flow rate (21) without the 
effects of moisture Separator. 


Initially assume state point 1 thermodynamic 
properties are ideal (21). 


Thermodynamic pump work (23). 


Tube profile, flow parameters across the tube 
bank (24, etc.). 


Tube sheet diameter (30). 


Evaporator shellside pressure drop for two phase 
flow (33). 


. Moisture separator pressure drop (538). 


Properties at state points 3 and 4 (39-41). 
Revised ammonia mass flow rate and velocity (50) 
includes the effects of the moisture separator; 
Mmecrate with (31), 


Revised ammonia heat transfer coefficient (5l, etc.); 
1terate with (9). 


. Heat exchanger cost analysis. 


In the following section, the basic steps summarized 


above will be described in detail. 


ANALYSIS OF THE EVAPORATOR AND MOISTURE SEPARATOR 


Salt Water Mass Flow rate, Msw 


The salt water mass flow rate through the hot pipe 


must be equivalent to the flow rate through the evaporator 


(assuming no leakage) 
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im, (Hor pie ) = iN (EvarorzATor) 


and mM = K3. AV eo 
where A = cross-sectional area of the hot pipe. 
V = salt water velocity through hot pipe. 
/B = density of salt water evaluated for an 


average hot pipe salt water temperature. 


As previously stated, the diameter of the hot 
pipe and salt water velocity are among the initializing 
conditions of the optimization and will be treated as design 
variables. 


2. Total Number of Evaporator Tubes, \} 


Using equation (1), it follows that 


: 2 
M, Si it ti V4 Ne (2) 
| 4 


where iia salt water density evaluated at the average 
bulk temperature initially assumed as the hot 


pipe salt water temperature, 


Zo 





oe = tube inner diameter. 
! : : : 
Ns = the number of tubes required to maintain 
the mass flow rate for an average salt water 


velocity per tube. 


The total number of tubes can be determined by solving 
mee (2) for Ne. 

The diameter of the tube and average salt water 
velocity per tube are initialized for the analysis and will 
be treated as design variables by the otpimization code. 


3. Total Evaporator Heat Transfer Surface Area (Outer), A; 


Having determined the number of evaporator tubes, the 
total heat transfer surface area can be determined using 
initializing values of outer tube diameter and tube length. 


For tubes without extended surfaces 


Ap = de Ly Ne (3) 


As previously, the outer tube diameter and tube length 
are initializing conditions and will be treated as design 
variables. 


4. Overall Heat Transfer Coefficient, [j 


The quantity "U'' provides a measure of the total 
thermal resistance in the flow path, based on either inside 
Or outside surface area. 

This analysis will be based on the value of U for the 


outside surface area derived from Eq. (3). 
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Using a resistance analysis, assuming one dimensional 


(radial) heat flow, 





Tsw OTe Tw1 Tw2 Ttnh, Ths 
_— © 
R 1 R 2 R Ra 


Rs 


the overall heat transfer coefficient may be expressed as 





Ae A. ee nclen «As (4) 
Th hw Ai "A iol ZK + Fenus + Ho hung 
where sw = tubeside heat transfer coefficient. 
Resw = salt water fouling heat transfer resistance. 
K = thermal conductivity of the tube material. 
Bel = outer and inner tube diameter. 


ate SIMO Oingeheab thanster resistance 
(assumed to be negligible). 
Wee outer and inner total fin efficiency (for 
plain tube analysis, total fin efficiency 
equals 1). 
A, = total outer surface area (including fin 


and bare tube). 
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A; = total inner surface area (including fin 


and bare tube).! 


where Agni; = total inner fin surface area. 
A; = total inner surface area (rien Gite im 
and bare tube). 
Agn, = total outer fin surface area. 
eee Oca Outer SUrtace area (Including fin 
and bare tube). 
Rs: = fin efficiency of single interanl fin. 
iz = fin efficiency of single external fin. 
a. Tubeside Reynolds Number, fre, 
Since the heat transfer coefficient correlations 
for the evaporator and condenser are dependent on tubeside 
flow, Reynolds number must be calculated. 


The tube Reynolds number is defined as 


Re, = EN SO (5) 


1Note that this analysis will hereafter consider smooth 
plain tube configurations only. 
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mere 2is,, = dynamic Viscosity of salt water. 
ene densteyeet Salt water. 


INtelaliyve properties are evaluated for 
Tauck = Tee (INLET) (6) 


Reynolds numbers greater than 2300 will be indica- 
tive of turbulent flow [Ref. 3]. Transition flow was considered 
laminar for numerical evaluation. 

Poste Natwemenedt: lranster, Coefficient, lane 

The simple empirical relation proposed by Sieder 

widelace [Ref. 3], expressed as 


5 3 0.14 


if | 
) ut 
Nuy= 18% (Fes Pe) (| ral (7) 
Le} \“4w 
was used for laminar heat transfer in tubes as defined by 
nem (5). 
Nusselt and Prandtl numbers, Nu, and f- , are 


defined as 
Nug = pects 
Kw 


i = Cai 430 
Ksw 


Where Lfsw > Cpy,, and ksi fdynamie: Viscosity, specific heat, 
and thermal conductivity) of salt water are evaluated at salt 


water bulk temperature. 


f 
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The Vemeecrsor the viscosity ratio term in 


Eq. (7) al 0.14 
vA 


where 4, is salt water viscosity evaluated at tube wall 
temperature, is considered negligible and will hereafter be 
dropped from the expression of Eq. (7). 
Relation (7) is based upon the following assumptions: 
fully developed flow in smooth tubes. 


fluid properties are evaluated at the 
bulk fluid temperature. 


and 1s valid for the following condition 
Ke, Pr A > 10 
E 
For fully developed turbulent flow in a tube as 
defined by Eq. (5), the Dittus-Boelter correlation [Ref. 3] 


expressed as 
3 
Nuy (5 Oa Re, a (8) 


was used. Nusselt and Prandtl numbers, Nuy and fr , are 
previously defined by Eq. (9). 
Relation (8) is based upon the following assump- 
Euons: 
fully developed flow in smooth tubes. 


fluid properties are evaluated at the bulk 
fluid temperature 


and is valid for the following conditions: 
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Prandtl numbers ranging from 0.6 to 100. 


. moderate temperature differences between 
Ciel leand tluid Conditions, 


c. Salt Water Fouling Heat Transfer Resistance 
In this document, it will be assumed that the 
fouling resistance coefficient for tubeside salt water can 
be maintained at .00025 (h-.t*F/aTu) using one of the 
following techniques: 
Chilerinat ion. 
MAN Brush System. 
Amertap. 
Chemical cleaning 
Pressure drops associated with cleaning techniques 
will not be considered in this analysis. Piping losses will 
be a function of tube length, inner diameter, salt water 
velocity and the absolute roughness of the tubing design 
material only. 
d. Ammonia Shellside Heat Transfer Coefficient, Niue 


Pek easally Nas will be assumed 


hung = 1000 (BTU/hr. ft" F*) (9) 


Since its value cannot be directly calculated during this 
phase of the analysis. 
ictneetne thermal Geolstance Expressed as 


Ri ee Ay 


Ay 


1} 





f 
Cro 
C2 <a 


a nee Ai 


\\ 


i 


ion ot ay 
Rs = ee 


2 


AG aes 


fieiatetial yalue for the overall heat transfer coefficient may 
bemcalculated. 
Ll. = : (10) 
en 
Rit R22 +RztRs 
5. NTU-effectiveness Relations 
The NTU-effectiveness relationships will be used to 
determine the evaporator outlet salt water temperature, 


Currently, all salt water properties have been based upon 


the initial assumption that 
Teux = Ty, (SW INLET To evar ) 


The expression for the number of transfer units (NTU) 
which is a measure of the size of the heat exchanger is given 


by 


NGG = legAe: / Gary, 


where Gaels defined as capacity rate of the single phase 


flow in an evaporative or condensing two phase flow regime, 
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Cimin = aac. Gre C2) 


Evaporator effectiveness can then be expressed as 


(-NTU) 
GE=f{-€ (Zz) 


for two phase flow regardless of the flow geometry. 
Using the definition of effectiveness 
actual heat transfer 


Effectiveness = —————————_______________ (ee) 
maximum possible heat transfer 


€ = cl _ Almin _ Tai - Tue (14) 
Omax f\T max Tees - Tc; 
The expression for /\7Tminrepresents the single phase 
(salt water) flow and Te: represents ammonia inlet temperature 
to evaporator taken at state point 3A. 


6. Evaporator Salt Water Outlet Temperature and 
u emperature 


Using the relationships of Eqs. (12) and (14), the 
following expression may be formulated for salt water outlet 


temperature 


(-NTU) 


Te = ie (Tui - Tez )(1- = ) (Se) 


GS: 





Concurrently, a revised evaporator average salt 


water temperature can be expressed as 


Tauck = ( Tt; i Tite )/2 (16) 


Using the revised value for average salt water 
temperature,iterate with equation (1) until the revised and 
current values of bulk temperature satisfy a specified 


convergence criterion. 


@ 


7. Amount of Heat Absorption, G 


Wsomewcne Gesults of Eq. (16) and (12), the amount of 


heat absorption by the evaporator may be expressed as 


on Cla = Tis) (17) 


8. Log Mean Temperature Difference, LMT D 
The NTU-effectiveness method can be used to determine 
the mean effective temperature difference (LM/)) across the 
evaporator (heat exchanger). 


Using Eq. (17) and the definition of 


G oe LI, A. E LMTD 
pec Pray = Gan A Tmax 


the log mean temperature difference across the evaporator 


may be expressed as 
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(-wTtu) 7 
( Ae a yer — fae) 


a (18) 
las Te 
where Te; = Jans evaluated at state point 3. 


i 


correction factor on LM7TD, equal to 1 for 
two phase flow. 


9. Film Temperature for Property Evaluation, [¢ 


In order to evaluate the shellside ammonia heat 
transfer coefficient, working fluid properties (i.e., 
viscosity, specific heat, etc.) must be evaluated at the film 
temperature to validate critical heat transfer expressions. 

By modifying the expression in Eq. (10) multiplying 
by a Single tube outer area, a value for single tube con- 


ductance can be expressed as 


A 


ul,A = = 
Sia ie ICS Ie aa as 


Subsequently, the average amount of heat transferred 


per tube would equate to 


O= LULA Gites, 


Led 


where T3 = Thats evaluated at state point 3. 


Again using the resistance analysis in Section 3, 


shellside wall temperature may be derived from 


Tae = Tuck = 2 a amd 
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Knowing Shellside wall temperature and the free-stream 
temperature, film temperature can be derived from their 


arithmetic mean. 


ie = ie ofa Ts 


(19) 
JE 


10. Ammonia Mass Flow Rate, nae 


According to first law of thermodynamics for steady 


State, steady-flow conditions in the evaporator: 


EVAPORATOR (ideal) Q(hot sw) 





Mens he +P = nas hs (20) 


from which the ammonia mass flow rate, Mas > may be determined 
if the enthalpies at state points 2 and 3 are known. 

If we initialize the lower and upper bounds of the 
analysis in terms of pressure Ps and Pe, respectively, and 
initially assume that a saturated vapor leaves the evaporator, 


the following relations may be expressed 
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hy = h,| P4 Th = Tsar} py 


oa) 
hs = hg) ps Ts = Tsar] P3 
where h4 = Tepresents enthalpy at state point 1 at 
EHewSUGEIomeIiniet to the working fluid 
Clee te Om PUldD . 
hs = represents enthalpy at (ideal)/state point 3 


as a Saturated vapor. 


Ty, 13 Seucuecce Mem eNenncopCetive Saturation 
Cempe BarUre on 
U4 = represents the specific volume at state 
peant 1. 


To summarize, the upper and lower pressure bounds of 
the system ( Pj and fs) will be initialized in the analysis 
and treated as design variables by the otpimization code. 
Temperature at state point 3 iS initially assumed to be a 
Saturated vapor (ideal T3); however, the working fluid is 
subject to a shellside pressure drop as it passes across 
the evaporator with an outlet quality of 90-95%. Properties 
at state point 3 (actual) will be assessed in follow-on 


sections. 


oF 





AMMONIA CIRC PUMP 





Meng D1 + Wer = Muys h2 (22) 


Assuming steady state, steady-incompressible flow, 
the change in kinetic and potential energies, and heat 
losses are negligible for isentropic conditions, and the 


isentropic pump work can be expressed as 
Wee. =) U4 Us oa FP) where Pz = fee 


After the isentropic pump work is calculated, the actual 


(adiabatic) pump work may be determined using pump efficiency, 


p: 





4} = WW Fee 
Wee = Wer. 
1p 


Actual outlet enthalpy at state point 2 may be deter- 


(23) 


mined uSing the results of Eq. (23) with Eq. (22) knowing the 
emenmaipy at state point 1 from Eq. (21). 

Heme the resultSsoL Eqs: (21) and (22), the mass 
flow rate in Eq. (20) may be calculated as the average shell- 


Side mass flow rate for the working fluid (ammonia). 
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11. Tube Profile, Flow across Tube Bank, and Tube Sheet 
Wizomte Ger 


Since the heat-exchanger arrangements (evaporator 
and condenser) involve multiple rows of tubes, the geometric 
arrangement of the tube profiles 1S important in the 
determination of the heat transfer coefficient, the tube 
Sheet diameter and the shell side pressure drop associated 
with two-phase flow (homogeneous model) [Ref. 4]. 


The following geometric arrangements are used: 


cone 
oe 


IN- LINE rs 
where Sn = pitch ratio x outer tube diameter, equal to Sp. 
pale = PutchieG wale. sencmiadls cance sbDetmech tube centers 


With respect to outer tube diameter. 


Ap = tube profile area (centerline to centerline) 
veh be. 
Se eee (24) 
pe 
Ap = Sn (25) 





Le Sp 


STAGGERED 30 Sy eg FLOW 


ta] 
2 FR ds sm sc (26) 


Sn 


Pp = en aaeeo > aie (27) 


Therefore, the tube profile area (centerline to centerline) 


per tube 1S equal to 
A» = Sn 29 (28) 


iesratio Of minimum flow area to the frontal area can be 


expressed as 
J Sn ay 


Using the selected tube profile geometry, either in- 
line or staggered, and knowing the required number of tubes 
by equation (2), the tube sheet diameter for heat exchanger 


design can be assessed as follows: 


40 





Zz 





Ni Ap = tl Tse 
ee 7 (30) 
where [so = tube sheet diameter. 

To estimate the shellside ammonia flow velocity the 
following control volume is introduced (ammonia circulation 
piping and the top portion of the evaporator). 

E 
i iis ini: a oar es eoae | 
ee | 
| 
| See Sera =! 
‘So? ceeereteeerer 
nat 

If the mass flow rate remains unchanged across any 

Peundary (continuity), 
ey Mee 

Furthermore, if we assume the evaporator has the 
means to evenly distribute liquid droplets across the top 
of the tube bundle (spray nozzles and baffling), the follow- 
i”emexpressions can be applied to estimate the mean droplet 
Mele@erty approaching the bundle: 

Let (Ag hig = Ag} 
where VY = percent of tube frontal area which is occupied 


by edrombe ts. 


The mass flow rates are 


M> = 4 Ap Vp 


4] 


ea 





Mg = (Aric Vi 


where Ap = ammonia pipe cross-sectional area. 
\/p = average ammonia velocity in the pipe. 
Therefore 
Ap 
Vp = A P 
(Ar) 
and since , 
eee 


At 
it follows that the average velocity of ammonia through the 


Circulation pipe is equivalent to the average velocity of 


ammonia at the tube frontal area boundary. 
Vo = V5 (31) 


Thus the assumption that K=Ap/As Ls MeN bve Werke ice 
the assumption of constant liquid kinetic energy in the 
Posmemwen rom the pipe exit to the bundle entrance. Con- 
Sidering the minimum free-flow area for shellside flow passage, 


Age can be derived from Eqs. (29) and (30): 


Bie = i L+ 


Age = Ag ae. cl) (52) 


De, 


= 
OR 
4) 
ry 
rq) 
css 
a 
| 


= represents the flow frontal area, 


aes 
ct 
I 


tube length. 
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Using the calculated values of Eqs. (32) and (20), 
the mass velocity for the minimum free-flow area can be 
expressed : 
ey are iNet 
where Myx represents the average ammonia mass flow rate. 


12. Pressure Drop of Two-Phase Flow across a Bank of 
Tubes, /\{ 


This portion of the analysis will use an analytical 
model for two-phase pressure drops applicable for a fog or 
spray flow pattern occurring at high void fractions -- the 
homogeneous model [Ref. 4]. 

The model asserts that if the pressure drop in the 
two-phase flow for a liquid-vapor mixture 1s relatively small 
compared to the absolute pressure, the flow is considered 
incompressible. Subsequently, the density of each phase is 
practically constant. During the process of phase change, 
the phase and velocity distributions are changed, and so is 
the momentum of the flow. Therefore, the pressure drop of 
a vertical two-phase flow consists of three components: 
friction loss, momentum change, and elevation pressure drop 
arising from the Pee of the gravitational force field. 

The local pressure gradient for a two-phase flow 


may be expressed as 


Be = Alpecront F yonenTUM + AFevevation (53) 


gene ig 
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For a given channel length, [L 


mie pressure arop 


@emponentesS Can be represented by 


and the total pressure drop, /\ 





dll 
| ag 
We - Sie a 
Ej | 34 
FRICTION Maun (34) 
2s 
Pm: 
CN igesiee een = Je 
AP. main <th 
ELEVATION ~— - C 
we 
P 
ee 1S given by the sum of 


these expressions 


where $ = 


Le = 


where X = 


Single-phase friction factor by Jakob expressed 
in Eqs. (35) and (36). 

GiriincimalOW tenet dertined for horizontal 
tubed evaporators as Le=lsp (tube sheet 
diameter). 

equivalent diameter of flow channel, defined 

by De = Fa dy - clo. 


mean specific volume defined by 


quality of mixture (State point 3). 
specific volume of liquid (state point 1). 


specific volume of vapor (state point 3). 


44 





The basic assumptions of the homogeneous model (fog 
flow model) [Ref. 4] are: 
(1) equal linear velocities of vapor and liquid, 
(2) thermodynamic equilibrium between the two 
phases, and 
(3) a suitably defined single-phase friction factor 
1s applicable to the two-phase flow. 
Using assumption (3) and the correlations by Jakob 
[Ref. 3], a suitable single-phase friction factor can be 
calculated from previously defined tube profile relationships: 


for staggered tube arrangements: 


-0-16 
ae OAL Ee 
f = O25 + Tc 7 3 =e Re. (55) 
Sin - ae / do | - 
and for in-line tube arrangements: 
Ost 
2a S | 
EY Ln oe (36) 
 ( Sn -de)/ao | CRRA 


Where Reynolds number (max) is determined from the shellside 
ammonia flow and the nozzling effect of the tube geometry as 


expressed by 


Mirae = V¢ s-) 


where Ve = the ammonia velocity at the tube frontal area 


boundary determined by equation (31). 
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Reynolds number for maximum shellside flow can be 


calculated using the following expression 


fe = Le Vmax Ble 


Ax Wace (37) 


Eq. (37) and tube profile data can then be used to 
evaluate the single-phase friction factor, required for 
Eq. (34). <All other components of the total pressure drop 
Eq. (33) can be determined from previously calculated data. 

13. Pressure Drop Across the Moisture Separator. APmsep 

This portion of the analysis will simulate the use 
of a cyclone separator to improve te evaporator outlet vapor 
quality. The flow pattern in a cyclone Separator is complex 
and simplifying assumptions are inadequate to allow the 
calculation of the corresponding pressure drop, which can 
Wmeryecrom | to 20 inlet velocity heads [Ref. 5]. Therefore, 
the worst case condition will be applied with an approxima- 
tion for the fluid flow inlet area to the separator banks. 

By approximating the inlet area as a fraction of 


the evaporator frontal area 


Aimer = Of Tso Ly 


the inlet fluid velocity can then be determined using the 


Wouking t1luid mass flow rate, Eq. (20). 
Dae = P Aimer V 


where Ve =e density of ammonia at State point 3. 
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Therefore, if the pressure drop across the moisture 


Separator is equal to 20 times the inlet velocity head, 





-~ 202 V 
OP sce wha ae (38) 


14, Enthalpy at State Points 3 and 4 
Since Eq. (33) represents the pressure drop across 
the evaporator shellside, the actual pressure at state 


point 3 or evaporator outlet may be determined from 


FS wew) s Ps -AP evap (39) 


where Fz was previously described as the pressure for a 


Saturated vapor. 


Similarly the actual pressure at state point 4, the 


moisture separator outlet, may be expressed as 
y= Pawew)—APy cep (40) 


Operating under the dome of the Temperature-Entropy 


diagram, the following properties are defined 


n3¢ (new) = h¢ | pseuew) hae = he] Pe 
(41) 


N3q (ew) = 1g | P3(vew) hag > hg) Pa 
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The subscript (NEW) representing a revised property 
Will hereafter be dropped from the expressions in Eq. (41). 
Assuming an evaporator outlet quality of 90-95%, and a 
moisture separator outlet quality of 99-99.5%, enthalpies 
at State points 3 and 4 may be determined using the relation- 


Ships of Eqs. (41) 
hz = har + X32 (h3q - hss) 


ha = has + X4 (haq- Las) (42) 


15. Revised Ammonia Mass Flow Rate and Velocity 
Till now, we assumed that the shellside mass flow 
rate was given in accordance with the ideal system defined 
by Eq. (20); however, in actuality this 1s not the case. 
The diagramatic representation that follows better 
illustrates the heat absorption phase of the OTEC power 
system and will provide the basis for the analysis and 
optimization. 
Note, aS in the previous control volume analysis, 
the following conditions are assumed. 
Steady state. 
Steady-incompressible flow. 


Change in potential and kinetic energies is 
negligible. 


48 













MOISTURE 
SEPARATOR 





EVAPORATOR 







GRAVITY DRAIN 


Analyzing the moisture separator as a separate 


control volume, 





If we assume that there is no carry-over of vapor in 


the separator drain, then 


and 


flow rate through the separator drain must be included in 


x35 Nz, = X4 Na 


M3 eee AAP 


However, for reasons of flow continuity, the mass 


the control volume analysis; therefore 


ng = M4 + Mp2 
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(43) 


C44) 





Substituting Eq. (43) into (44) and solving for /N% e-, 


the following expression can be derived 


a) A So 
Me = (“x3 |) ma (45) 


Looking at the evaporator as a separate control 


volume, 





m@emeneTesy balance is 


Gt tizh2 + Mriher + Ht: hor 2 0133 + Mee bre 


Assuming the change in enthalpy across the re-flux 
pump and the difference between the separator drain outlet 
and evaporator inlet are negligible, the energy balance 


becomes 


Grme2h2 + Moe Noe = 1113 Ng (46) 
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where hoe= he || paces) fluid drained from the separator is 


assumed to be a saturated liquid, 
Furthermore, a mass balance of the evaporator control 


volume can be expressed as 
M2 + MR. te Mo; = M3 + Mee (47) 


Where WkKe = MR; 


s0iving Eq. (47) for the mass flow rate at state 
point 3 and substituting into Eq. (46) with Eq. (45) yields 
the following expression 
5. m2h2+/x4 . = ye cae hs 
In addition, a mass balance for steady-state, steady-flow 


indicates that the mass flow rates at state points 2 and 4 


are equal and therefore 


a= M2 (49) 


Using Eqs. (48) and (49), the revised mass flow 
rate at state point 2 may be determined. Concurrently, the 
revised average ammonia velocity acting on the tube profile 
geometry may be determined from this revised mass flow rate. 

Using the revised ammonia velocity acting on the 


tube profile geometry and iterating from Eq. (31) until 
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an acceptable convergence criterion is achieved provides 
the pressure drops across the evaporator and moisture 
separator, and the properties at state points 3 and 4 for 
a given film temperature. The result is more representative 
of the heat absorption phase in the OTEC power cycle than is 
the commonly used ideal analysis. 

In addition, solving for the revised temperature at 


Seace POlnt 3, 


Tz = Tsar | ps (es) 


and iterating through Eq. (18) revises the film temperature 
and subsequent working fluid properties. 
16. Revised Shellside Ammonia Heat Transfer Coefficient 

In the search for acceptable correlations to predict 
the average evaporative heat transfer coefficient, two analyti- 
cal treatments were found that lent themselves to OTEC 
power system conditions. 

The first of these correlations seeks to predict 
Elimetiim eVaporation heat transfer coefficient for horizontal 
tubes [Ref. 6]. Owens [Ref. 6] uses (1) the similarity 
between evaporation and condensation, (2) the correlation 
forms of local evaporation heat transfer coefficients for 
water on a vertical tube developed by Chun and Seban, and 
(3) the dependence of heat transfer on the vertical spacing 
of the tubes as was experimentally demonstrated by Liu, to 
arrive at the following correlations for non-boiling thin 


film evaporation: 
i 512 





for laminar flow 


Q.1 WES =—t/5 


h = 22 (4 Ce (46 | (51) 


Gouetcurbulent flow 


= Cue . 7 Ovo 
h= tes (1) at (cect (52) 
a S/S lage K¢ 


where H = vertical Spacing with respect to tube outer 





diameter. 
i. 


The laminar-turbulent transition point is defined 


tube flow rate per unit length. 


by the intersection of Eqs. (51) and (52) 


is 


~ 1490 Ce 4+ 
He, = 168 ( a) 


K F 


The pseudo-Reynolds number for horizontal vertical falling 


film evaporation is defined by Ref. 7. 


Re= 40 
Alf 


The second correlation combines boiling and evapora- 
muon or liquid films om horizontal tubes, applicable for 


vertical banks of plain and enhanced tubes [Ref. 8]. 


oe 





The overall model for a single tube is expressed as 


—_ 


h= hy+ha be + he ({- fa | (53) 
L 


where aye = Rohsenow pool boiling correlation over the 


entire tube length given by 


Can y= VFF es 
at 1 Csp = function of the fluid-surface combination. 


Z\l = wall temperature minus free stream saturation 


temperature. 

Wee = Surface tension 

ha Sloe ehanes FoumcOc ttt Clentermmtne developing 
region. 





mA4/S ; 
Ly = | ( earn ; 
41s \ ary 
and Ae = fully developed heat transfer coefficient given for 


laminar flow by 


| sae -1/f5 70.22 
he = 082 (- ae (55) 
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and, for turbulent flow, 


e175 O.4 SOF 








-3 ia ie og ae 
he = 38X10 (3 (se ye | ($6) 
a er 
where L, = circumferential length of heated surface. 


thermal diffusivity. 


Lg 


[" = flow rate per unit axial length of tube. 


developing length around tube circumference. 


Hemappiy Eq; (51)})) for a vertical bank of tubes, L 


is expressed as 


i N, ido /2 


The laminar-turbulent transition point is defined by 


the intersection of Eqs. (55) and (56) 
— {.06 


eo Pe 
ee - 5800(= 


As before, the pseudo-Reynolds number is defined by Ref. 7 


ee, oh le (57) 


After using Eq. (57) to establish which flow regime 
the system is operating in, the revised heat transfer coeffi- 


cient for non-boiling thin film evaporation or nucleate 
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boiling may be calculated and then iterated with the initial 
assumption for the shellside heat transfer coefficient, Eq. (9). 
This will have a convergence effect on variables which are 

a function of the shellside heat transfer coefficient, moving 
them closer to actual OTEC system performance characteristics. 

The user should be aware that the predictions for the 
OTEC power system using ammonia have been for the case where 
no boiling occurs in the film. This condition is dictated 
by industrial preference for plain tube heat exchangers to 
Minimize fouling and the characteristic of ammonia to wet 
surfaces well, flooding out nucleation sites. A number of 
enhancement techniques have been developed to create 
nucleate boiling, including a variety of tube configurations 
and surface preparations; however, a preference for them 
has not materialized. The nucleate boiling development in 
Eq. (51) which would be indicative of tube enhancement is 
provided for information only and will not be included in 
the optimization or summary of conclusions. 

Havangedesecribed the Methods used to predict the 
Shellside heat transfer coefficient, we can complete this 
chapter of the OTEC power system analysis by constructing the 
heat exchanger cost analysis. 

17. Evaporator Cost Analysis 

At the request of TRW, Wyatt Industries, a large 
exchanger fabricator, prepared cost estimates for three 
different sizes of vertically configured evaporators and 


condensers, based upon initial design specifications prepared 
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by TRW. Based upon these estimates, TRW developed sets of 
equations that represent the costs of various heat exchanger 
component parts for shell diameters ranging from 10-35 ft 
ama 55-50 ft [Ref. 9]. 

The following are the TRW evaporator cost ($) 
equations as a function of outer tube diameter (inch), 
total number of tubes and tube-sheet diameter (ft) for 
tube-sheet diameters of 10-35 ft. 


Drilling time/tube sheet thickness 


tas 0.66 (ad, - 0.5) (58) 


. Thickness of the tube sheet 


0.68 
= 0.56 Is, (59) 


. Tube sheet labor cost 


z 156695 (Nz /9630)( tg /0,06) te, /4) (60) 


Orsi 


. Tube sheet material cost 


— == 


pote 2.5 


. Tube installation cost 


0.7 


. Heat exchanger shell cost 


Ghyess 177265 (Lyt6\ Tro /ie) (63) 


Si 





. Ammonia distribution plate and battles cost 
7 2 ae yf } / . | = 
C npg = 19865.79 (Ne/ae30\(, / 0.66 }( Tsp /18) (64) 
Bustle, flanges channels and flow plates cost 
2 
Capee = 3OBSSO (T/18) (65) 


. Tube material cost 


Cry = (5 L. +z) Nz Ae (66) 
tS 
where C1 = curve fit of tube cost per foot. 
£2 = tube machining cost if required 
. Heat exchanger head costs 
= 
Co ay 253240 (Ts0/18) (67) 


. Water inlet, nozzles and supports cost 


| 2 
Cums = 220510. 15 (Tsp /18) (68) 
Tube welding costs (Titanium tubes) 
for Ne < 26000 
ele) O.7 
Cry = 1473 Ne (ado/1-5) ae 


For N, > 36000 | By; 
f 13 _ 
The sum of cost Eqs. (60) through (69) would equal 


the cost to fabricate one OTEC evaporator with a tube sheet 
diameter of 10-35 feet (all the preceding component costs 
have been adjusted for current pricing at a 10% annual rate 


Sie tua f10n ) . 
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If our analysis is based on a 30-year life-cycle 
criterion, no adjustments are necessary to any component 
cost equation if titanium tubing 1s used due to its anti- 
corrosive qualities; however, using aluminum tubing (l.e., 
A1l-5052), the expense of retubing must be considered to meet 
the criterion of a 10-year life cycle for aluminum tubing. 
This implies Eq. (61) and (66) must be modified to reflect 


the costs of retubing at the 10 and 20-year point in the cycle. 


. Aluminum tube installation cost 
SIS , 20 
Care = Crr[te ei) + Ue | (70) 


where l = projected inflationary rate (input by customer) 


. Aluminum tube material cost 
20 
. \ {0 ; 


For tube sheet diameters of 35-50 ft the following cost 


relationships apply [Ref. 9]: 
Equations for drilling time/tube sheet 
thickness (58), thickness of tube sheet (59), 
and tube material costs remain unchanged. 


Tube sheet labor and material cost (titanium) 


0.141 0:68 

Coe = ODEN Gy tae (72) 
2,014 

Cism 7 29.566 Isp Tel (73) 


£5) 





Tube sheet labor and material cost (aluminum) 


Odi 0.68 
? 


{,61 
ee = BS, Tso Cr. es) 


. Tube installation costs 


Ca 
Cp = 56.542Ni ado (76) 


. Heat exchanger shell cost 


2-06 
Give = {2.544 (Le +6) Tsp (77) 


. Ammonia distribution plate and baffle costs 


1.82 0.873 
Cp, = 158.079 Isp +7Z2.419N, tt, (78) 


. Bustle, flanges, channels, and flow plate costs 


— Bie 
Coes 2 EV TYNES (79) 
. Heat exchanger head cost 
LAS 
Cryy = 1125-3! thes (80) 


~ Water inlet, nozzies and support cost 


— 


_ cee fe 
Crene = (CS aa Isp (81) 
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Tube welding costs (titanium tubes) 


Belg Ne & NOS, 


mos ec 
Cry = 14.73 Ne (do /15) (82) 


for Ny > 36000 


(03> N, 
Cay = ONE Ne ao/us)” 


As indicated previously, the cost to fabricate one 
OTEC evaporator with a tube sheet diameter 35 to 50 ft is 
equal to the sum of component costs Eqs. (72) through (83) 
(all the preceding component costs have been adjusted for 
current pricing at a 10% annual rate of inflation). 


For an analysis based on a 30-year system life-cycle 





criterion, the additional costs for aluminum retubing must 


be considered and Eqs. (70) and (71) apply. 
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bees PARASITIC LOSSES 


A. INTRODUCTION 
This chapter describes in detail the programming analysis 
for parasitic losses which include: (1) pumping and pipe 
requirements for both cold and hot salt water systems, 
(2) pumping and pipe requirements for the working fluid 
(ammonia) circulation and re-flux systems, and (3) turbine 
generator losses due to inefficiencies. Hotel requirements 
have not been incorporated into the analysis, but could be 
included for the final design analysis. 
Pumping power requirements will be determined through 
the use of the general energy equation between the inlet and 
outlet of the system control volume [Ref. 3]. 
2 2 : 
AE VE gg Bi = Me gg Bet We + (bovses) 
Va Ze ras 


4 


To determine the pumping power Wg, the following effects 
will be evaluated: 
1. Density head. 
Zeer rictionmbhosses. 
Intake piping. 
. Heat exchanger tubing. 
Exit piping (if employed). 
3. Thermodynamic pressure head. 


4. Elevation head. 
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See tiner Losses. 
Intake piping inlet configuration (contraction). 
Intake piping Screen (obstruction). 
PlOWeetnoupievalves, elbows, etc. 
Outlet palpinow(expansion). 
Inlet to heat exchanger tubing (contraction). 

. Qutlet from heat exchanger tubing (expansion). 

Outlet of exit piping (if employed). 

In the above pump head evaluations, the following inputs 

are specified: 


Pipe lengths (hot, cold, ammonia circulation and 
re-flux piping). 


iimiem pape diameters (initialized and treated as a 
design variable by the optimization code). 


. Absolute roughness corresponding to piping/tubing 
material (designer specified). 


Fluid velocity (initialized and treated as a 
design variable by the optimization code). 


Pump mechanical and electrical efficiencies. 

As an overview of the parasitic pump loss analysis, the 
following major steps in the algorithm are listed in order 
of their execution: 

. Hot pipe salt water pump. 
Inlet piping friction losses (86). 


. Minor piping losses due to inlet Screen (87) 
and plenum design to evaporator core (88). 


Evaporator core minor losses (89, 90) and 
tubeside friction losses. 


Total pressure losses (92) and pumping 
head (93). 


Pumping power requirements (95). 
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Pump cost analysis (96). 
Cold spawewsatrc Water pump. 


Initialize cold pipe inner diameter and SW 
velocity (design variables). 


. Minor losses due to inlet ducting (97) and 
plenum design to condenser core (98). 


(Tike teplpingunmterlon losses (99). 


Condenser core minor losses (100, 101) and 
tubeside friction losses (103). 


Density head (104). 
Total pressure losses (105). 
Pumping power requirement (107). 


Pump cost analysis (108). 


. Ammonia circulation pump. 


. Piping friction (109) and minor losses due to 
valving/elbows (110). 


Pressure drop across evaporator shellside (112). 
Thermodynamic head (113). 

Elevation head (114). 

Total pressure losses (115). 

Pumping power requirement (116). 


Pump cost analysis (118). 


. Ammonia re-flux pump. 


Piping friction (119) and minor losses due to 
valving/elbows (120). 


Thermodynamic head due to pressure drop of 
Saturated liquid ammonia across evaporator 
Shellside (122). 

Elevation head (123). 


Total pressure losses (124). 


64 





.. Pumping power requirements (126). 
Pumomeost amalysis (127), 
Parasitic pump losses. 
In the following section, the basic steps summarized above 


Will be described in detail. 


PeeeoaNALYSIS OF PARASITIC LOSSES 
i Hot Pipe Salt Water Pump, up 
The pressure losses due to piping friction and 
associated minor losses will be determined using the 


Darcy-Weisbach correlation [Ref. 10]. 


Nn ; a . 
A\P= > papas cen 
teh Ze 


where Kj; describes the resistance coefficient. 


<< 
r 


fluid velocity. 


: L 
Ki = f= (84) 


where mimic ClLONataccor. 


equivalent length in pipe diameters. 


£ 
D 
In order to determine the friction factor, the pipe 


flow Reynolds number must be calculated. 


2 Josvi Vsn Ai 
Ke, Faas Sw 
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where ey Sa properties of walt water at the hot pipe 
inlet temperature (assumed constant through- 
Mitmene Digme) . 
acne salt water velocity and inner pipe diameter 
(initialized and treated as design variables 
by the optimization code); velocity assumed 
constant over pipe length. 
Pipe flow Reynolds number greater than 2300 will be 
considered turbulent. 


for laminar flow 
f. . 
a 85 
ey 
for turbulent flow 


poraen 


Re ery (86) 
[in (E/5.7 di + 5.74/Rey 4 |” 


where € = absolute roughness corresponding to piping material 
selected. 
Eq. (86) yields a friction factor within one percent 
of the Colebrook equation and is valid for the following 


conditions [Ref. 9]. 
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Considering the resistance coefficient for pipe minor 
BoSSes 


. Assume the inlet duct is the same size as the 
pipe inner diameter, but it is screened 


Ke ee (87) 


. ASsume piping enters evaporator through an 
area which is abruptly changed [Ref. 11] 


< [ a (ait) ] (88) 


where ley = evaporator tube sheet diameter (assume tube 
Sheet diameter is twice as large as the inner 
pipe diameter). 

SUmIMineg themnesules OLM@eas. (S4)e (87), and (88)mto 
determine the total resistance coefficient, the pressure 
losses due to piping can then be determined using Eq. (83). 

If a variety of valves or fittings are to be included 
with Eq. (84), Ref. 11 provides a representative listing of 
equivalent length-to-pipe-diameter values. 

To analyze the pressure drop across the evaporator 
tubeside, we again use the Darcy-Weisbach correlation, but 
for different design assumptions. 


Assume inlets to evaporator tubing are well 
rounded [Ref. 11] 


kz 0.5 —) 
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. Assume outlets of evaporator tubing expand 
to an infinite reservoir [Ref. 10] 


Ke (90) 


Using the Reynolds number in the previous chapter, 


meeei 5), the corresponding friction factor Eq. (85) or (86), 


and resistance coefficient can be determined 


where 


io EQS . 


Bn 
log 2 = evaporator tube length and inner tube 


diameter and are initialized and treated as 
design variables by the optimization code. 
Summing the results of the resistance coefficient 


(89), (90) and (91), the pressure losses due to the 


evaporator design may be determined using the Darcy-Weisbach 


@ermreiation Eq. (83). 


The results of the piping losses and core design 


losses are equivalent to the hot pipe salt water pumping 


system requirements 


A Prime = APowve sysren + APevap vesien (92) 


converting to pumping head 


Ge | 
H = Bag Ot pame (93) 
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Pumping power in terms of horsepower can be determined 


using the following expression 





new | gH 
fs = Msvw ( q (94) 
MG : 
where Ne = pump mechanical efficiency (designer input). 


Msyj * salt water mass flow rate determined in 
previous chapter, Eq. (2). 
To equate parasitic pump losses to power input, 
Eq. (94) is converted to the motor load requirement in terms 


of megawatts electrical. 


FAP mw) = Pue x CONVERSION FACTOR (95) 
M 


where TAG = pump motor efficiency (designer input). 

Because of the high salt water flow rates and rela- 
tively low pumping heads, good engineering design would 
dictate the use of axial flow (propeller) type pumps. 

Using the algorithm developed by TRW [Ref. 9] from 
data provided by Johnston Pump Co., and Process Equipment 
Co. (distributors of Ingersoll Rank and Johnston Pumps), 
the cost of salt water pumps can be expressed as 


2 
C (D/1000) 0.75 + Se et Xa (96) 


pump ~ 
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where 


d= Zz ve, 


where A; Vsy = inner hot pipe diameter, salt water velocity 
(initialized for analysis and treated as 
design variables by the optimization code). 
The above algorithm is valid for the following con- 
ditions 
vertical, wet pit, propeller type pumps with 
cast lron steel columns with protective epoxy 
coating, stainless steel shaft and bronze 


impeller. 


. pump size from 155,000 through 750,000 GPM 
With total dynamic heads of 8 through 12 feet. 


Eq. (96) has been adjusted for current pricing at a 10% 
annual rate of inflation. 
ge COld Pipe Salt Water Pump, Ke 


Using Reynolds number 


R: _ fsa Vigne 
ed ~ LS 


where (sw fe, * properties of salt water at the cold pipe 





inlet temperature (assumed constant through- 
out the pipe). 

Vsw Ai = salt water velocity and inner pipe diameter 
(initialized and treated as design variable 
by the optimization code), velocity assumed 


constant over pipe length. 
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Pace tlow characteristics and friction factor can be identi- 
fied. A pumping analysis will be developed for the cold 
pipe pump using the Darcy-Weisbach correlation, similar to 
the development in the preceding section. 

Considering the resistance coefficient for minor pipe 
mOsSes 


. Assume the inlet duct is well rounded [Ref. 11]. 


K= 09 (97) 


iN@e lf 


. Assume piping enters condenser through an 
area which is abruptly changed [Ref. 10]. 


=| 1+ (di a (98) 


ee mm 


where ‘lap = condenser tube sheet diameter (assume tube sheet 
diameter 1S twice as large as the inner pipe 


diameter). 


. Assume one ninety-degree elbow is required 
in system [Ref. 11]. 


L 30 
D 
Summing the results of Eqs. (84), (97), and (98), 


the total resistance coefficient can be expressed as 


ie L , 
K= 5 ( vn r 5) . K ike i K pseu vee 
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where L =wlenethmorecold pipe. 

A; = inner diameter of cold pipe. 

Pressure losses due to piping can then be determined 
uSing the Darcy-Weisbach, Eq. (83). 

In analyzing the pressure drop across the condenser 
tubeside, the Darcy-Weisbach correlation is used again, but 
for different design assumptions. 


. Assume inlets to evaporator tubing are well 
rounded. 


k=0.5 


(100) 
. Assume outlet of condenser tubing expands to 
an infinite reservoir. 
K=10 (101) 
Defining Reynolds number for condenser tubeside 
flow, while assuming 
Tauck = | Cono Cincer) (102) 


= w Vsw d: 
i 


esi 


where Jars , so = properties evaluated at condenser tubeside 
bulk temperature (initially assumed equal to 


cold pipe inlet temperature). 
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Vsw, 4i = average salt water velocity through 
tubing, inner condenser tube diameter (both 
are initialized and treated as design 
variables by the optimization code). 
mie cOrrespOnding friction factor, Eq. (85) or (86), and 


resistance coefficient can be determined 


Naya = if Le (103) 


Zhe 


where Ai, the condenser tube length and inner tube diameter 
are initialized and treated as design variables by the optimi- 
gaeton code. | 

Summing the results of the resistance coefficient 
in Eqs. (100), (101), and (103), the pressure losses due to 
the condenser design may be determined using the Darcy-Weisbach 
@eprelation, Eq. (83). 

A complete analysis of cold pipe losses must also 
include the effect of density head and a corresponding increase 
in pumping power requirements. 

For most engineering problems involving the flow of 
liquids through a pipe, where the temperature change in the 
pipe is small, the density of the fluid is considered to be 
a constant and the fluid is termed "incompressible." However, 
the flow problem in OTEC cold pipe systems is unique. We can 
continue to assume that there is negligible change in the fluid 
temperature, virtually unaffected by the ocean thermal 


gradients, because of the system's characteristic high mass 


iO 





flow rates. However, the height of the water column (1500 
to 3000 feet) inside the pipe requires the effect of fluid 
compressibility to be taken into consideration. 
The effect of an increase in density with depth 

can be expressed by the integral 

€ 5 
mnpoceid 
Vas) 

d 

With a density head defined as? 


£ a 
Hp SE rf = 


a 





Integrating the pressure-density variation, the 


density head reduces to [Ref. 12] 
Ho= 22-Fi + A_ (P.-P.)|1- Kon (Pe + P:) | 
Z fg a 


where Kn = mean compressibility of walt water, f(salinity, 
temperature and pressure). 
reference density at which K,, is evaluated. 


f 


Considering pressure at any depth obtained from the 


mmcegral , 


P= ~ gf 242 


the density head can be rewritten as follows 


Be fe 


Ze 
=! BIT: opal as — Km =) 
He = (2c Z ) af poe | Ke [ og 


*Note that Z is measured as positive upward so that ocean 
depth values (Z., 2: ) are negative and G25) 1s a positive 
quantity. 


| 
| 


74 





Rigorous procedures for calculating the density 
procile which is a function of temperature, salinity and 
pressure may be found in Ref. 13; however, they will not 
De discussed in this document. 

For the purposes of simplification, the following 
solution technique was developed: 

(1) If the liquid in the pipe is taken to have a 
constant density with respect to pressure, the compressibility 


approaches zero; the density head can then be expressed as 


He = (222) 5 [70 


(2) Converting the geometric term for elevation to 


an equivalent integral expression 


Ze 


a ee | ie 
ai eran zl 74 


The reference density is taken to be the inlet value so that 


it 


and the density head can be rewritten as follows 


Zz 
se (2: - 2(2)) a2 


(3) Assuming a linear distribution of density with 


depth, due to temperature variations, as illustrated below 


ies 





—————E——E—————— ee coeeeee”rcTeeeeeee_eerOoOvwvre--- 





Ze wap =a! 


the following linear expression for density with respect to 


depth may be formulated, where 2,=© for convenience. 


fia > ain~pe)U- 2/2i) 


(4) Applying the equation developed in section 3 to 
the density head integral above and integrating over the 
range of values for sea water depth (z), the following 
equation is derived as a linear approximation to the density 


variation of sea water with respect to depth 


= (Pi. = fre ; 
7 (2 -) 
where (Pi = curve fit evaluations of density for 
specified depths of sea water. Data 
extracted from Ref. 14, 
The results of the piping losses, core design losses, 
and density head are equivalent to the cold pipe salt water 


pumping system requirements 


A Foump = Df pire System + LF cong DESIGN + AP pens iry (105) 
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emcee quecomeeca. (105) can be converted to a 
pumping head. Similarly, pumping power in terms of horse- 


power can be determined using Eq. (44). 


een as (2 
= 


where 





Ms = Ves (Ha Ves ee) 

(ONG 

and josn = density of salt water evaluated for a 
constant inlet temperature. 

Vsw Ai = cold pipe salt water velocity, and inner 
diameter (initialized and treated as design 
variables by the optimization code). Note 
Salt water velocity through cold pipe is 
considered to be constant. 

Pumping power can then be expressed in terms of 


megawatts electrical 


feed) = Pep x CONVERSION FACTOR exe) 
M 
where “ves pump motor efficiency (designer input). 
Using the same arguments for the selection of an 
axial flow (impeller type) pump, as used for the hot pipe 
salt water pump, the pump cost algorithm developed by TRW 


can be applied to the cold pipe salt water pump assuming 
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the required conditions are validated. 


Cum = 


( D/1000) 0.78 +80] 1.21 x10" ae 
Equation (108) has been adjusted for current pricing at a 
10% annual rate of inflation. 
3. Ammonia Circulation Pump, ire 
The function of the ammonia circulation pump is to 
circulate and lift saturated liquid ammonia from the condenser 
hot well at state point 1 and increase its pressure to exceed 
the operating conditions in the evaporator at state point 2. 
In order to evaluate these characteristics, the 
following pumping elements will be included in the analysis: 
Piping losses (friction and minor). 
Heat exchanger shellside pressure drop. 
Thermodynamic pressure head. 
Elevation head. 
As in the preceding analysis, Reynolds number is used 
to determine pipe flow characteristics 
Reg= Ay 


Saturated liquid properties of ammonia for 


where a 


the temperature at state point 1 (assume any 
temperature increase from pump work is negligible). 
ai = inner pipe diameter (initialized and treated as 


a design variable by the optimization code). 
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\/ = ammonia flow velocity determined from the pre- 
ceding chapter, Eq. (50). 
The ammonia pipe friction factor can then be deter- 
mined from Eqs. (85) or (86), and the piping friction 


resistance coefficient can be expressed as 


i = fe (109) 


where | = ammonia circulation pipe length (designer input). 
couse rine ties nesistcance Coerficient for miner pipe 


losses, assume there are four ninety-degree elbows in the 


system 


Cnet) 


where = 3 equivalent length in pipe diameters for a standard 
elbow [Ref. 11]. 
Summing the results of Eqs. (109) and (110), piping 
losses (friction and minor) can be determined using the Darcy- 


Weisbach equation (83). 


* ic 2 
AFape= PL f4 +44 “= a 


The heat exchanger shellside pressure drop is also 
included in the pumping head requirement because it serves as 


a resistance to flow. 
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Pressure drop across the evaporator shellside was 
determined using the two-phase flow model (homogeneous) 


expressed by Eq. (33) 


A ap = Aer en +L ouentun + Ol erevarion (112) 


Since the pump is required to lift the working fluid 
to a higher elevation and increase its operating pressure, 
the following elements must be included in the analysis: 


Thermodynamic head 


where = a a Jiggs 

L anne Fez ( ) 
represents the difference in thermodynamic 
operating pressure between state point 2 and 
State point l. 


Elevation head 


INST A Frrevarion= 72771 (114) 


ZA 
2a 


datum. 


il 


elevation of the evaporator inlet above 
datum (taken to be equal to evaporator 
tube sheet diameter plus 25). 

represents the lift head required to move the 


working fluid to a higher elevation. 
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The results of piping losses (111), evaporator 
pressure drop (112), the thermodynamic head (113) and 
elevation head (114) are equivalent to the ammonia circulation 


pump system requireements. 


AFoump = L\ Pore te DN ee NE AA aero +A\ FeLevaTion (115) 


Using Eq. (93) with ammonia properties, Eq. (115) 
Can be converted to pumping head and finally expressed as 


pumping power (horsepower). 





get ( | (116) 


where lege = mass flow rate of ammonia determined by 
Eq. (20) of the previous chapter. 
Nop = pump mechanical efficiency (designer input). 
Pumping power can then be expressed in terms of mega- 


Watts electrical 


Ferre x CONVERSION FACTOR 7) 
ie 


fene(mw) = 


where ae pump motor efficiency (designer input). 
Because of high pumping head and moderate flow rates, 
good engineering design would dictate the use of a single 


Suction centrifugal flow type pump. 
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Using the algorithm developed by Westinghouse Electric 
Co. [Ref. 15] from data provided by Bingham Pump Division, 
Portland, Oregon, the cost of the ammonia circulation pump 


can be expressed as 


; _ ed 5 
G meninges | 7p IC 
purnp — | “75” fx NC (118) 
BOltOo 
where Ns = mass flow rate of ammonia (1b, /he) 
Vg = specific volume of saturated liquid ammonia 


at state point l (Ft*/Ibn) 
Eq. (118) has been adjusted for current pricing at a 10% 
annual rate of inflation. 
4, Ammonia Re-flux Pump, Pae-riux 

The function of the re-flux pump is to recycle 
ammonia droplets which are not evaporated in the heat absorp- 
tion process. Saturated liquid at approximately the heat 
exchanger's operating pressure is lifted from the evaporator 
drain to the ammonia feed inlet, for redistribution as droplets 
across the evaporator tube bundle. (Drainage mass flow rate 
1s assumed to be equal to 30% of the evaporator inlet feed 
mass flow rate.) 

In order to evaluate these characteristics, the 
following pump elements will be analyzed: 


Piping losses (friction and minor). 


82 





Thermodynamic pressure head. 
Elevation head. 
As in the preceding analysis, Reynolds number is used 
to determine pipe flow characteristics 
Re, - = 


saturated liquid properties of ammonia for the 


where Fie 
average pressure across the evaporator. 
di = inner pipe diameter (initialized and treated 
as a design variable by the optimization code). 
VY = ammonia flow velocity determined from the 
evaporator drainage mass flow rate assumed 
equal to 30% of the evaporator inlet feed mass 
flow rate (assume velocity constant throughout 
the pipe). 
The re-flux pipe friction factor can be determined 
from Eqs. (85) or (86), and the piping resistance coefficient 


can be expressed as 


Ke fh (119) 


where L = ammonia re-flux pipe length (designer input). 
Siecedsilnee considering the resistance coefficient 
for minor pipe losses assume there are four ninety-degree 


elbows in the system 
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(120) 


where L = equivalent length in pipe diameters from a standard 
elbow. 
Summing the results of Eqs. (119) and (120), piping 
losses (friction and minor) can be determined using the 


Darcy-Weisbach, equation (83) 


L yl 

cl f£i4gh ( s 12a 

AF; 08 yal ae iu D ZW ( 
In order to determine the thermodynamic pressure 

head, the pressure drop across the evaporator for the 

Saturated ammonia liquid must be analyzed. Since the 

Saturated vapor and liquid are in thermodynamic equilibrium, 


the results of Eq. (112) apply. Therefore 
AF. = Ps -P2 


Therefore, the thermodynamic pressure head is equal 
to the pressure drop across the evaporator for the saturated 


ammonia liquid. 
AP ryeamo = Sug (122) 


Finally, the elevation head is equal to the elevation 
of the evaporator feed inlet with respect to datum, the 


@rain outlet. 
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Therefore, 


A Fey 2 OLS) 
where Zoe dacune drain outlet. 
Z2 = elevation of the evaporator inlet above datum 


(taken to be equal to the evaporator tube 
Sheet diameter plus 10). 
The results of piping losses (121), the thermodynamic 
pressure head (122), and elevation head (123) are equivalent 


to the ammonia re-flux pump system requirements. 


AF 


; ao le (124) 


THERMO 


PiPe ELEVATION 


As before, using Eq. (93), Eq. (124) can be converted 
to a pump head and finally expressed in terms of pumping 


power (horsepower). 





eure (3H) (125) 
Mp \ Je 
where Ma = drainage mass flow rate. 
No = pump mechanical efficiency (designer input). 


Pumping power can be expressed in terms of megawatts 


electrical 





FRE FLUK yy) FREAK, CoyyERSION FACTOR (126) 
M 
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where Neo = PUMpmMobowectiicaency (deSigner input). 
Using the same arguments for the selection of a 
centrifugal pump, the pump cost algorithm developed by 


Westinghouse can also be applied to the ammonia re-flux 


pump. 
: 0.04 e 
- {Max Ve) 421 

(i - yaglex 10 

NE (Fae Glan) 
where Me = mass flow rate of evaporator drainage ammonia 

(1b,,/hr) 

Vz, = specific volume evaluated at the average 


evaporator pressure (ft ji | bm} 
Eq. (127) has been adjusted for current pricing at 
a 10% annual rate of inflation. 
5. Parasitic Pump Losses 
Parasitic pump losses is the summation of electrical 
auxiliary pumping requirements (hotel and maintenance loads 


not included) determined by Eqs. (95), (107), and (126). 


Foss > Pup i Pep ue Vemee r Poewerux (128) 
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Voom Ne AND ELECTRICAL POWER 


A. INTRODUCTION 

The turbine generator is one of the critical elements 
of the OTEC power system. Its energy conversion efficiency 
and efficiency of design have a major effect on the overall 
System performance. To illustrate this point, Ref. 16 
reported that a three-point change in turbine efficiency from 
85 to 88% results in a 3.6% increase in gross power, and a 5% 
increase in net power developed. 

This chapter will describe the analysis to evaluate the 
expansion turbine thermodynamic properties and generator 
output. The use of these properties will determine the 
internal turbine efficiency and outlet quality subject to 
design and thermodynamic constraints. The relationship 
between the condenser operating pressure (design variable) 
and the turbine outlet quality will be used to initialize 
the heat rejection characteristics of the condenser. 

General literature on turbomachinery designed for OTEC 
closed cycle systems indicates that a turbine having the 
following characteristics 

Double flow, axial inflow, 
Four stages of expansion, 
Operating at 1800 RPM, 
provides the optimum aerodynamic design [Ref. 16]. However, 


it is not the intent of this thesis to analyze the geometry 


Sy 
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and performance parameters of the turbine. Turbine geometry 
such as 
Specific speed and specific diameter, 
eiiccmciecdmne tek. 
Rotational speed, 
Blade height, 
Blade stresses, 
Should be treated as a separate systems problem using 
optimization to improve state-of-the-art design. 
Parasitic losses due to the following generator turbine 
inefficiencies will be evaluated in this section. 
Generator mechanical and electrical. 
Turbine mechanical. 
As an overview of the turbine-generator analysis, the 
following major steps of the algorithm are listed in order 
of their execution: 


Gross electrical output with no parasitic losses 
leZ9)... 


Enthalpy at state point 5 (130). 
Turbine outlet quality (131). 


Entropy at state point from a specified outlet 
Celi ese TOL 4 a 


Quality amd entha@ipy at state point Ss» (133, 134). 
Internal (adiabatic) turbine efficiency (135). 
Turbine cost analysis (137). 
Generator cost analysis (138). 

In the following section, the basic steps summarized above 


Will be described in detail. 
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MeeeANALYSIS OF THE TURBINE AND ELECTRICAL POWER REQUIREMENTS 
mee GEOss Electrical Output and Inefficiency Losses 
[If the net electrical output required is indicated 
by (in terms of megawatts), the gross electrical load at 


the turbine shaft can be expressed as 


= ie 
Ver Mgen 
where Pee = parasitic pump losses determined by Eq. (128). 
oes) cUgoine Mechanical efficiency (designer 
input). 
Neen = generator mechanical and electrical 


efficiency (designer input). 
The loss of electrical output due to generator- 


turbine inefficiencies is equal to 


o eo { 
OSs V4 Neen 
2. Turbine Efficiency 


The power developed across the turbine is 
Fo=m™ (hs-ha) 


where m mass flow rate of ammonia given by Eq. (48). 


h4 


From this, the enthalpy at state point 5 can be calculated. 


enthalpy at state point 4, Eq. (42). 


If we initialize the operating pressure of the con- 


denser in terms of Ps, the following relations may be expressed 
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hsq = hg )ps hse = he | Ps (130) 


Therefore, it follows that the turbine outlet quality, 


X5, can be determined from 
hs = hss xs (hsq- hse) (isi) 


Having established the moisture separator outlet 
pressure and temperature, Eqs. (40) and (41), the entropy 
at state point 4 can be determined for a known separator 
outlet quality (designer input) using the following 


relations 


Saf = Se] T4 S4g = Sq) fea 


54 = Sap + X4 (Say a S4¢) (132) 


For isentropic turbine work, 
54-2558 (133) 


the quality at state point 5s may be determined using the 


following relations 
59 = S4)rs Ss¢= Se) nee 
Ses 3 Sst X55(Ssq zs Ses) (134) 


Having determined the quality at state point 5s, the 


enthalpy can now be determined. 
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hss = hss + X55 (hsq- hse) Csi 


De tnigettemnescmumeamor Eqs. (41), (130), and (132), 
the internal turbine efficiency (adiabatic) can be deter- 


mined, expressed by 


h4 - hs 


Vee ea (136) 


ha - hess 


litmeiolincwaerealistic selection of internal efficiency, 
the following constraints are attached to the optimization 
code 
» hs < hs 
Cu OO 


3. Turbine Cost Analysis 
The ammonia turbine cost is based on an algorithm 
developed by Westinghouse to estimate manufacturing costs 


eet. 15). 


bf | oe KG 
a = 2.42X10 (0375+ Ee /136000 N;) Fe (137) 


GEoOssvetectrical Output in KW. 


where ey 


2 (for a double flow turbine). 


Pad 
ae 
u 


Fe = flow price factor (1.0 for single-flow, 1.447 
for double-flow). 


oA. 





The above algorithm is valid for the following 
conditions: 
Double flow, axial inflow. 
PS VMawlti-stage. 
Operating at 1800 RPM. 
The generator cost will be based on an algorithm 


developed by TRW from data provided by selected manufacturers, 


| : é 
Coew = (0.023 E, +063) 1,21x10 (138) 


and is valid for the following conditions 
1800 RPM rotor speed. 
. power factor 0.8. 
Eqs. (137) and (138) have been adjusted for 


current pricing at a 10% annual rate of inflation. 
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VI. CONDENSER 


A. INTRODUCTION 

As indicated in the introduction to Chapter III, several 
heat exchanger concepts have been proposed for the closed- 
cycle OTEC system, with variations in their design. 

The analysis to be presented for the condensing heat 
exchanger will be based upon the following design 
characteristics: 

Single-pass shell and tube heat exchanger. 


Horizontal/vertical orientation of tubes with an 
equilateral triangle or square tube profile. 


Smooth plain-tube configuration (no enhancements). 


Tube material (titanium or aluminum based on a 
30-year life-cycle criterion). 


Biofouling control based upon an achievable 
OUMMANO ide GO 


Heat exchanger centerline located on sea surface. 
As an overview of the condenser analysis, the following 
Major steps in the algorithm are listed in order of their 
execution: 
Initialization of design variables (DV). 
Tube length. 
SW velocity through condenser tubes. 
. Outer tube diameter. 
Tube profile pitch ratio. 
Amount of heat rejection (139). 


Tubeside bulk temperature (142). 
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miotel mumber Of tubes (143). 
Log mean temperature difference (144). 
Conductance (146). 
Number of transfer units (145). 

. Heat exchanger effectiveness (147). 


Initially assume a value for ammonia heat 
transfer coefficient (151). 


Single tube conductance (148). 
. Average heat rejection per tube (152). 
Film temperature (153). 


Revised ammonia heat transfer coefficient (154, etc.); 
iterate with (151). 


Tube profile, flow parameters across the tube 
bank (158, etc.). 


. Tube sheet diameter (163). 


Condenser shellside pressure drop for two-phase 
flow (166). 


ReviscdepmhOperties at State point 1 (171, 172); 
iterate with (21). 


Overall heat transfer coefficient (173). 
Total heat transfer surface area (174). 
Revised condenser tube length (175). 
. Heat exchanger cost analysis. 
In the following section, the basic steps summarized 


above will be described in detail. 


B, ANALYSIS OF THE CONDENSER 


i mount of deat Rejection, @’ 


Using the calculated value for enthalpy at state 


point 5, equation (131) from the previous chapter, the ideal 
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values at state point 1, Eq. (21), and the steady-state mass 
flow rate of ammonia, Eq. (48), the amount of heat rejected 


by the condenser can be expressed as 
G = Mins (he-hi) (139) 


2. Tubeside Bulk Temperature 
As in condenser tubeside Reynolds number, salt 
water properties will be evaluated at bulk temperature, 
initially assumed equal to the cold pipe inlet temperature. 
Using this premise, the condenser salt water capacity 


rate can be evaluated 
Cimin = (Nae CD. (140) 


where Cae specific heat of salt water initially evaluated 
at the cold pipe inlet temperature. 
Mep = mass flow rate of salt water through the cold 
pipe previously evaluated by Eq. (107). 
Using the results of Eqs. (139) and (140), and the 
known cold pipe inlet temperature, the condenser salt water 


Outlet temperature may be evaluated from the basic expression 
G = Cimin (Te,- ieay (141) 


where Jc, Tc, = condenser salt water outlet and inlet 
temperatures, respectively. 
Having determined the condenser salt water outlet 


temperature, the revised bulk temperature can be expressed as 


a5 





Ta = Fee lev (142) 
A 

Using the revised condenser bulk temperature and 
iterating with Eq. (102) corrects the operating temperature 
for salt water properties which are essential to the 
analysis. 

3. Total Number of Condenser Tubes, Ny 

Since the mass flow rate of salt water through the 

cold pipe is equivalent to the mass flow rate through the 


condenser, according to the law of continuity, 
Mep = Memo 


it follows that the number of condenser tubes for a specified 


tube diameter, can be evaluated using the following expression: 


n= ou La V4 N, (143) 


where [sw = average salt water density evaluated at 


bulk temperature. 


d; = inner tube diameter (initialized and treated 
as a design variable by the optimization code). 
Vi = average salt water velocity through the con- 


denser (initialized and treated as a design 


variable by the optimization code). 
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4. Log Mean Temperature Difference, IMD 
Using the result of Eq. (141), the known pipe salt 
water inlet temperature, and the inlet temperature of 
ammonia evaluated at state point 5, the LM/D of the condenser 


may be expressed as 
iMi p= es 1ci (144) 


5.  NTU-Effectiveness Relations 
The number of transfer units which 1s a measure of 
the condenser size can be determined from the basic 
expression 
Nru = We Ae (145) 
Cmin 
where the conductance (ULAs) of the heat exchanger is a 


Munetion of the heat absorbed and the LM7TD. 


The condenser effectiveness can then be expressed as 


(-NTU) 
€= {-& G47 | 


for a two-phase flow, regardless of the flow geometry. 








6. Single-Tube Conductance, L.A. 
Using the resistance analysis derived in Chapter III, 
Section 4 for an initialized tube length 
Lee 
the heat exchanger conductance for a single tube can be 
expressed as 


HAs = . 





{ | je (148) 
— k.. In d./« a t in PST 
Mi Nsw Ai Ai foi ZirKL Ag nity t ie Ro 


tubeside heat transfer coefficient. 


where nee 
Kesw = salt water fouling heat transfer resistance. 
K = thermal conductivity of the tube material. 
A.A; = Oral Olucermaid singer tube Surface areas 
(including fin and bare tube); tube length 
1S initialized and treated as a design 
variable by the optimization code). 
Rents = ammonia fouling heat transfer resistance 
No Ni = outer and inner total fin efficiency 
a. Tubeside Reynolds Number 


Since the salt water heat transfer correlation 


is dependent on tubeside flow, Reynolds number must be 


evaluated R. - ya Via 
=a ae 
fis 
where [30/450 = salt water density and viscosity are 


evaluated for the fluid's bulk temperature. 


| 
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oh Vea = inner diameter and average salt water tube 
velocity. 
Reynolds numbers greater than 2300 will be 
indicative of turbulent flow [Ref. ais 
Decale Waeer Heat Transfer Coefficient, ew 
Once again the empirical relationship proposed 
by Sieder and Tate [Ref. 3] will be used for laminar heat 


transfer in tubes and as defined by 


Ve VE 


Nluj= 1.86 (Re, P) (4 (a 


Ww 


Nusselt and Prandtl numbers are defined as 


Nu y= Rew di (149) 


P= CPs Mw (150) 
Ksw 
where dynamic viscosity, specific heat, and thermal conduc- 
tivity of salt water are evaluated at the salt water bulk 
temperature. 
The effect of the viscosity ratio in the 
Sieder-Tate equation is considered negligible, and will 
hereafter be dropped from the expression. The assumptions 
and validity condition associated with the Sieder-Tate 
equation were stated in Chapter III, Section 4, and will not 


be repeated here. 


sls 





roOueetwbyeaaveloped turbulent flow, again the 


Dittus=Boelter correlation [Ref. 3] was used. 


Nua = 0.025 Pf; = 

Nusselt and Prandtl numbers are previously 
defined by Eqs. (149) and (150). Assumptions and conditions 
for validity were stated in Chapter III, Section 4. 

c. Salt Water Fouling Heat Transfer Resistance 

As indicated previously, it will be 
assumed that the foulding resistance for tubeside salt 
water can be maintained at 00025 (hr. f¢-F/atu) 


d. Ammonia Shellside Heat Transfer Coefficient, hyy, 


inital eleva Nes will be assumed 


Dine = 1600 (BTU/hr.#t: F) (151) 


Since its value cannot be directly calculated during this 
phase of the analysis. 
Using the following single-tube thermal 


resistance j 


R,= ————__—___ 
Mj Ae dil 








LOG 





an initial value for single tube conductance (outer tube 


surface) may be calculated 


{ 
L.A. = —-——--——_ 
R,+Re2+ Ra+ RS 


7. Film Temperature for Property Evaluation, [¢ 
In order to evaluate the shellside ammonia heat 
transfer coefficient, working fluid properties must be 
evaluated at the film temperature. 
This can be accomplished by using the results of the 
Single tube conductance, the tube side bulk temperature and 
the working fluid saturation temperature, expressed in the 


following equation for single tube heat transfer rate (average). 


Q-= le (Ts Tauck) Gla) 


Again using the resistance analysis as in Chapter II1l, 


the shellside wall temperature may be expressed as 
Twz = |eutk +O(Ry + R24 Rs) 


Knowing the shellside wall temperature and the free- 
stream temperature, the film temperature can be derived from 
their arithmetic mean 


Te = Two+ Ie (153) 
ee 


For purposes of this calculation, saturated tempera- 


Mire COnditlions at state point 5 are taken to represent 
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free-stream conditions, when in fact the two-phase process 
will experience a pressure drop and a corresponding drop in 
temperature. 

8. Revised Shellside Ammonia Heat Transfer Coefficient, 


ut 


This analysis will include correlations for both 
horizontal and vertical heat exchangers. 

In the horizontal-tubed condenser, Nusselt's 
correlation was used as a predictor [Refs. 7 and 17], 


for laminar flow 


_ — {/3 
h- 045 ( ae, (154) 
yl Us 
where Ww = estimate of ammonia mass flow rate across 
each tube. 
Kee LU = properties evaluated at film temperature. 


L tube length (initialized and treated as 

a deSign variable by the optimization code). 
This correlation is probably conservative, since it 

does not consider turbulence due to high vapor velocity 

or splashing of condensate [Ref. 7]. 


For turbulent flow, Nusselt's correlation is increased 


by 10% as recommended by Jakob [Ref. 17] 


h = 1.045 ( K ; fe 4) (155) 
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The laminar-turbulent transition point is defined by 
a Reynolds number of 2100, where the pseudo-Reynolds number 
for film-type condensation on horizontal tubes is defined 


as [Ref. 7] — 
oa 
Vane 


where | = mass flow rate of condensate per tube over its 
length. 
In the vertical tubed condenser, both Nusselt's and 
Kirkbride's correlations were used as Deeductors ~| Retr. 7 |. 
For laminar flow, Nusselt's correlation is increased 


by a factor of 1.28 as recommended by McAdams [Ref. 7]: 





, } -1/3 -1/3 
i Jelg Aa 
ead a e ia teas 
y= 1.28 11.4 a (ao (156) 


where [= mass flow rate of condensate per tube over its 
diameter. 


For turbulent flow, Kirkbride's correlation is applied 


: 4 
= pe 
byw 0071 (| aie 
pe ovo eas) (Fe) ‘eal 


The laminar-turbulent transition point is defined by 
a Reynolds number of 1800, where the pseudo-Reynolds number 
for film-type condensation on vertical tubes is defined as 


[Ref, 7] 
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After using the pseudo-Reynolds number to establish 
the flow in which regime the system is operating, the revised 
heat transfer coefficient for film-type condensation may be 
calculated and then iterated with the initial assumption 
for the shellside heat transfer coefficient, Eq. (151). 

Once again this will have a convergence effect on variables 
in which the shellside heat transfer coefficient is a 
function, moving closer to actual OTEC system operating 
point characteristics. 


9. Tube Profile, Flow across Tube Bank, and Tube 
eet Diameter 


Since the condenser tube bundle involves multiple 
rows of tubes, the geometry of the tube profile arrangement 
1s important to determine the shellside heat transfer coeffi- 
cient, the tube sheet diameter and the shellside pressure 


drop associated with the "homogenous" two-phase flow model 


[Ref. 4]. 
Using the same arrangements shown in Chapter III, 
Section 2, Sp 
n 


Sn = Fa ae , (158) 


A pz on (159) 
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where Sn pitch ratio x outer tube diameter. 


rR 


pitch ratio (initialized and treated as a 


design variable for the optimization code). 


Ap 


tube profile area per tube 


Le Sp 


STAGGERED 30 S, —<g FLOW _ 


where 


Sym 23 2 IER ele rasta 


(160) 
Ap = 2n Sp (162) 
the ratio of minimum flow area to the frontal area can be 
expressed as 
Ase Seale 
a rs (265 
As aT; 


Using the selected tube profile geometry and know- 
ing the number of condenser tubes by Eq. (143), the tube Sheet 
diameter for the condenser design can be evaluated from the 


following expression 





| 2 
N. A» a tits (164) 


where Tsp = Tube sheet diameter. 
To analyze the shellside ammonia flow velocity, the 


following control volume is introduced (turbine generator 


discharge and top portion of the condenser). 


a 


| - 


! | 


“ee ip ss righ hhh had : 


Since the mass flow rate remains unchanged across 


any boundary 


% 


Furthermore, if we assume the condenser has the 
capability to evenly distribute vapor across the tube bundle 


(distribution baffles), the following development applies 


to the vapor coverage: 


fet (As) = AcY 


where K = percent of tube frontal area which is covered 


by vapor. 
Me = -% A. Vs 
M; = 5 Ag Ve i 
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where As Seeconcenser inlet cross-sectional area. 


Vea turbine discharge ammonia velocity. 


Therefore Ve = 45 As Vs 


(eAEy 


lf YsGAs/A » it follows that the turbine discharge 
velocity is equal to the average velocity of ammonia at the 
tube frontal area boundary. A determination of the distri- 
bution fraction / requires a detailed knowledge of the 
design of the turbine/condenser interface. In the absence 
of this information it is assumed that 

Ve = We 

A similar argument could be presented for a vertical 
tubed condenser where turbine discharge is admitted to a 
distribution ring that bands the condenser tube bank. 
Exhaust vapor would travel radially through the tube bundle 


megethen collect at the bottom after vertical film- 


i a 2 
condensation. — | 
—— 
al | 
l 
| | 
Vside aa 5 
— : 
a 
a ) 
~<—I ' 
i ae 


Again, in the absence of a detailed design, it 1s assumed 


feniiant 


Veeee Vs 
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Considering the minimum free-flow area for a horizontal 
tubed condenser, Ac¢ can be derived using Eq. (163) and the 


projected frontal area. 
Ag=lsp Le 


Ay=A, (-Zacde) aes 


the flow frontal area. 


= 
> 
4) 
sa 
4) 
a= 
ae 
i i! 


tube length. 
For vertical condensers 
Ac = ll Tso * FronraL LeneérH OF VAPOR INLET Flow 
Using the previously calculated value of the ammonia 
flow rate and Eq. (165), mass velocity for the minimum free 
flow area can be expressed as 
ae ths 


Aft 





(166) 


10. Pressure Drop of Two-Phase Flow across a Bank of 
uoes, AP 


The pressure drop in the two-phase flow condensing 
heat exchanger will be determined using the homogeneous 
model introduced in Chapter III. The model will consist of 
three components -- friction loss, momentum change, and 
elevation pressure drop arising from the effects of gravity. 

The local pressure drop for a two-phase flow may be 


expressed as 
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ae 7 L\ FRICTION ution (167) 


For a given channel length, ae the pressure drop 


components can be expressed by 


where 


—_ 
_ FtQ@aGy 
A FRICTION PDe2 Fe Le (168) 


co- 


ENP = at (169) 








MOMENTUM - Sle 
an 
Been U Ja = (170) 


Single phase friction factor by Jakob expressed 
amare (55 atom (56). 

mass flow velocity determined from Eq. (166). 
channel flow length, defined for horizontal 
tubed condensers as Lee Iso (tube sheet 
diameter) and for vertical tubed condensers 

as ee Lee (tube length). 

equivalent diameter of flow channel, defined by 


De = FR d,—- Ao 


mean specific volume defined by 


ae x 
aay [ 5 ee eC! )] 
= Ui { j- 
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where aaa Ooh Mixture (State point 5). 

feomepective volume of liquid (state point 1) 
oye t2esvolume Ct Vapor (state point 5). 
All components of the pressure drop model Eqs. (168, 169, 

and 170) can be determined using the preceding information. 

11. Revised Properties at State Point 1 

Since Eq. (167) represents the pressure drop across 

the condenser shellside, the actual pressure at state point l 


or condenser outlet may be determined from 


Fiwew)= Fi- APeouy (171) 


where py 1s previously described as the condenser operating 
pressure for the ideal cycle. 

Operating on the saturated liquid line on the 
Temperature-Entropy diagram, the following properties 


are defined: 


Se 
hy (NEw) = he | Pt (new) Tt (NEW) = Tsar] P4 (NEw) (GZ) 


The subscript (NEw) representing a revised property 
Will hereafter be dropped from the expression in Eq. (172). 
Until now, we assumed the condenser outlet tempera- 
ture and pressure were designed to operate as an ideal system, 
without a pressure drop. Therefore, using the revised 
temperature at state point 1 and iterating over the range 
from Eq. (21) until an acceptable convergence criterion is 


achieved, all the preceding variables as function of qh 
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Will be reevaluated to complete the closed-loop cycle of the 
Simulated OTEC power system. 
m2. Overall Heat Transfer Coefficient, | |, 

The quantity "U'"' represents a measure of the total 
thermal resistances in the flow path. Therefore, using the 
tube conductance expressed in Eq. (148) which is divided by 
the outer heat transfer surface area of a single tube, the 
overall heat transfer coefficient for the condenser can be 
determined. 


The thermal resistances are now expressed as 


_ A. 
i a ea “a0 


ils 
R 2 a ae 
Mi Aesw ae 


Rez = dolndo/ai 
IK 


Re= —~— 
Ve Awe 


and the overall heat transfer coefficient for the condenser 






may be calculated using 


, 
a 
a Ss tS (173) 


‘] 


lil 








moe ToegalGondenser Heat Transfer Surface Area, A: 
Having determined the corrected number of condenser 
transfer units (145), salt water capacity rate (140) and 
overall heat transfer rate (173), the total condenser heat 


transfer area can be calculated from the NTU expression 


Cmin ot) 


14. Revised Condenser Tube Length 
Using the total heat transfer surface area calcu- 
lated from Eq. (174) and the total number of condenser 
tubes (143), the revised condenser tube length can be deter- 


mined from the basic expression 
A.= INE Wale L; (Reviseo ) CS?) 


At this time, it 1S necessary ot iterate the condenser 
design until the two values (initial and revised) of the 
tube length converge. This iteration may be accomplished by 


the COPES routine if the following constraint is defined 
L pice ie Ly (Revisep)- Le (twtr at) 


Minimization of this difference will cause continual 
adjustment of the required tube length, already treated as 


a design variable by the optimization code. 
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ioe COndenser Heat Exchanger Cost Analysis 

As indicated in Chapter III, TRW developed sets of 
equations to represent the costs of various heat exchanger 
component parts for shall diameters ranging from 10-35 feet 
and 35-50 feet [Ref. 9]. 

The following are the TRW component cost equations 
for the condensing heat exchanger. Prior equation reference 
numbers will be substituted where equalities exist with the 
evaporative heat exchanger component cost expressions. 


for tube sheet diameter 10-35 feet 


Drilling time/tube sheet thickness. (58) 
Thickness of the tube sheet. (59) 
Tube sheet labor cost. (60) 
Tube sheet material cost. (61) 
Tube installation cost. (62) 
. Heat exchanger drill cost. (63) 


. Ammonia distribution plate and baffles 
cost. 


2 -2 Ze 

Cyog = 1539x10 ty Ne Top (176) 
Bustle, flanges, channels and flow 
Dilue seOst. 

2 20 
= S, ne 

ome {12 6 Isp Cee 
Tube material cost. 

ory = (C4 Ly oD) Na Bee (178) 


eS 





— 
=~ 





where C{ = curve fit of tube material cost 
per foot. 
C2= tube machining cost if required. 
Heat exchanger header cost. (67) 


Reid betmeli lc mmm Oiazese alld SuUppemt cost. 
Ginn = 10106.475 thes (179) 


Tube welding costs (Titanium tubes). (69) 

The sum of the preceding costs would equal the cost 
to fabricate one OTEC condenser with a tube sheet diameter 
of 10-35 feet (all the preceding component costs have been 
adjusted for current pricing at a 10% annual rate of 
mie2ation). 

If' our analysis is based on a 30-year life-cycle 
criterion, no adjustments are necessary to any component cost 
equation if titanium tubing is selected. However, using 
Al 5052-0, the expense of retubing must be considered to meet 
the 30-year life-cycle criterion, as in the cast of the 
evaporation. For convenience, and possible subsequent 
modification, these considerations are repeated here. 

Based upon the utility of Al 5052-0, two complete 
condenser retubings will be required to meet the basic 30-year 
criterion. This implies Eqs. (62) and (178) must be modified 
to reflect the costs of retubing at the 10 and 20 year point 


in the cycle. 





. Aluminum tube installation cost. 


10 20 
Care = Crr| 1+ (tri) + (trl) | 


where | = projected annual inflationary 


orem inpuee by Gls tomer ) . 


. Aluminum tube material cost. 


10 Rs 
Cor = Cam (tte) e(tri) | 


for tube sheet diameter 35-50 feet. 


Drilling time/tube sheet thickness 


. Thickness of the tube sheet. 


. Tube sheet labor and material costs 


(titanium). 


. Tube sheet labor and material costs 


(aluminum). 


. Tube installation cost. 
. tube material cost. 
. Heat exchanger shell cost. 


. Ammonia distribution plate and 


baffles cost. 


0.973 
Copp 7 9625 Ne te 


Bustle, flanges, channels and 
flow plate. 


2.184 


Heat exchange head cost. 


- 145 
= - 6z Y 
Cre 437 Tsp 


iS 


(180) 


(181) 


(58) 
(S9) 


(e257 5) 


(Ore 5) 
(76) 
(178) 
(77) 


(182) 


rs3) 


(184) 





Water inlet, nozzles, and supporters 


GO Sitee 
1.056 
Come = 1453.6 Is, (185) 
Tube welding cost (titanium tubes). (eo) 


As indicated previously, the cost to fabricate one 
OTEC condenser with a tube sheet diameter of 35 to 50 feet 
1s equal to the sum of component costs (note, all the 
preceding component costs have been adjusted for current 
pricing at a 10% annual inflation rate). 

For an analysis based on a 30-year life-cycle 
Criterion, the additional costs for replacing aluminum 


tubing must be considered and Eqs. (180) and (181) apply. 
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VII. NUMERICAL OPTIMIZATION 


A. INTRODUCTION 

Nearly all design processes attempt the minimization or 
maximization of some parameter or design objective. For 
the design to be acceptable, it must satisfy a set of con- 
straints which impose limits or bounds on design parameters. 

For the stated problem a computer program can be written 
to perform the basic analysis of the proposed design. If 
any parameters fall outside the prescribed bounds, the design 
engineer changes the parameters and re-runs the program. In 
effect, the computer code provides the analysis with the 
engineering making the actual design decisions. 

A logical extension to the computer-aided approach is a 
fully automated design, where the computer also makes the 
actual design decisions and performs trade-off studies. The 
COPES program provides this automated design and trade-off 
capability by the use of the optimization program COPES/CONMIN 
Heer. 18]. COPES is an acronym for Control Program for 
Engineering Synthesis, and CONMIN is an acronym for CONstrained 
function MINimization. Subsequently, a FORTRAN analysis 
program simulating a closed-cycle OTEC power system can be 
coupled to the COPES program for automated design, using 


_ some basic programming guidelines [Ref. 18]. 
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B,. COPES/CONMIN 

There are many numerical optimization schemes available 
to the engineer. Methods employed by these schemes fall 
into four basic categories: random search, sequential 
unconstrained Minimization, optimality criteria, and direct 
constrained optimization. The optimization progran, 
Selected for automated design analysis of the simulated 
OTEC power system, is based upon direct constrained 
optimization. 

Before any discussion of the optimization technique, 
basic definitions are summarized for convenient reference 
[Ref. 19]: 

Design variables - those parameters which the 
Optimization program is per- 
mitted to change in order to 
improve the program. 

Objective function - the parameter which is to be 
Minimized or maximized during 
Che wOpteImazatlon process. 

Inequality constraint - one-sided conditions which 
must be satisfied for an 
acceptable design. 

Equality constraint - condition which must be 
equaled for the design to be 


acceptable. 
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Side constraints - upper and lower bounds in a 
design variable, 

Assuming that the FORTRAN analysis program has been 
developed and a particular objective function has been 
Selected, the general optimization problem can be stated 
oe [|Ref. 20]: 


Find the vector of design variables,YX , to 


fen ze ely, (186) 


Subject to the constraints: 


G(x)<0 J=1,NCON (187) 

H:(x) =O J=1,NEQ (188) 

VLBSX:SVUB,  L=1,NDV (189) 

where X = the vector containing the set of independent 


design variables. 
Zege- Maicmoujectiue  tunceloneto be Minimized. 
G(X) = inequality constraint (NCON is the number 
Cee stem CONS tbaints) . 
Hj (x) = equality constraint (NE&@ is the number 
of such constraints). 
VLB;/VUB.= lower and upper bounds, respectively, 
On the design variables. 
If all inequalities of Eqs. (187) and (189) are satisfied, 
the design is said to be feasible if any constraint is not 


Satisfied, the design is infeasible. If the objective function 








is a minimum and the design is feasible, it is said to be 
the optimal design. 

meorder tO Staresthe @ptimization algorithm, the 
initial set of design variables, X, must be specified. I[t 
is desirable, but not essential, that the initial design 
Variables provide a feasible solution. The optimization 
algorithm will then proceed in an iterative fashion using 
the following relationship 


a a 


-— rd 


where ¢ = the iteration number. 
Xx = scalar quantity which defines the move in the 


search direction. 


way 
Ht 


vector search direction which will reduce 
the objective function (useable direction) 
Without violating constraints (feasible 
direction). 

To solve this problem, the optimization program 
COPES/CONMIN is used [Ref. 18]. CONMIN uses the Fletcher- 
Reeves algorithm for locally unconstrained problems [Ref. 20] 
and Zoutendijk's method of feasible directions (modified to 
improve efficiency and reliability and to deal with designs 
which do not initially satisfy all the constraints) for 
locally constrained problems [Ref. 21]. 

However, CONMIN does not handle equality constraints 
aorectly, but WEG by means of penalty parameters. To 


Pemreve thts. the .Opjective function 1s augmented as follows 





—— 


meet. 19]: 


Neg 


F (x) =i OUR KS H (190) 
ie 


and the equality condition of Eq. (188) is treated as an in- 
equality constraint 

H;(X) SO y=, NEP 
The penalty function approach effectively satisfies the 
equality constraint while maintaining the rapid convergence 
Characteristics of the CONMIN program. 

The numerical optimization problems of equations (186) 
through (190) are very general, allowing for any number of 
design variables and constraints. In assessing the value of 
optimization, the automated design provides a very attractive 
Seproach tO numerical optimization; however, there are both 
advantages and limitations to these techniques [Ref. 20]. 

Advantages: 

Reduction in design time. 
Systematic design procedure. 


. Applicable to a wide variety of design variables 
and constraints. 


Virtually always ylelds some design improvement. 

Not biased by engineering experience. 

Requires a Minimal amount of man-machine interface. 
Limitations: 

Computer times may increase dramatically as the 

number of design variables increases. A practical 


limit imposed by the current state of the art for 
most problems is 30 design variables. 
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Optimization techniques have no stored experience 
to draw upon; the validity of the result is 
limited to the validity of the analysis program, 


The results of the optimization are as correct 
as the analysis program is theoretically precise. 


Optimization algorithms used here cannot deal with 
discontinuous functions. 


The optimization program will not always obtain a 
global design optimum and may require restarting 
from several different points to acquire reason- 
able assurance of obtaining the global optimun, 
The analysis program must be properly structured 
to couple with the COPES/CONMIN optimization code. 
C. DESIGNATED DESIGN VARIABLES, CONSTRAINTS AND OBJECTIVE 
FUNCTION 
To assist in the interpretation of the enclosed OTEC 
power system FORTRAN analysis, the following summary identi- 
fies the design variables, constraint functions and objective 
function used in the analysis and subsequently operated 
upon by the COPES/CONMIN optimization code. These parameters 
are all contained in a labeled COMMON block in the computer 
code, referred to here as "GLOBAL COMMON.' Specific 
GLOBAL COMMON location numbers and upper/lower bounds for 
operating parameters summarized below can be located in 
Appendix C. 
Design Variables 
Inner cold pipe diameter 
Inner hot pipe diameter 
Inner ammonia circ pipe diameter 


Inner ammonia re-flux pipe diameter 
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Evaporator operating pressure 
Condenser operating pressure 

Outer condenser tube diameter 

Outer evaporator tube diameter 
Evaporator tube length 

Condenser tube length 

Condenser tube salt water velocity 
Cold pipe salt water velocity 
Evaporator tube salt water velocity 
Hot pipe salt water velocity 
Evaporator tube profile pitch ratio 


Condenser tube profile pitch ratio 


Goons traintePunctions 
Operating system pressure ratio 
Upper temperature bound of ammonia 
Lower temperature bound of ammonia 
Satisfactory enthalpy at state point 5 
Satisfactory quality at state point 5 
Satisfactory condenser tube length 
Internal turbine efficiency 
Evaporator tube sheet diameter 


Condenser tube sheet diameter 


Objective Function 


Cost of major power system components 
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VIII. CONCLUSIONS AND RECOMMENDATIONS 


ewe CONCLUSIONS 

1. The use of an analysis code for OTEC power systems 
coupled to COPES/CONMIN optimization code provides a power- 
ful tool to design an optimum power system for the desired 
net electrical output, measured against the objective 
mmmetion. Such a design could permit construction of higher 
capacity systems using the optimized modules as substations 
of the total power plant. 

2. The analysis code coupled to COPES/CONMIN provides 
an excellent vehicle to evaluate proposed designs relative to 
a true optimum. Tables 1 through 4 illustrate the result of 
preliminary calculations using the analysis code with an 
objective function to Minimize system cost. From these, the 
following conclusions can be drawn concerning horizontally 
oriented aluminum (A1-5052) and titanium-tubed heat exchanger 
power systems: 

a. The cost/KW output is nearly constant over the 
range of optimum designs for both titanium and aluminum tube 
heat exchangers. 

b. During testing for feasible plant designs in 
increments of 5 MW (net) electrical output, it was observed 
that a higher megawatt output plant could be achieved with 
titanium-tubed heat exchangers than for aluminum (A1-5052). 


For titanium-tubed heat exchangers, a 25 MW (net) power 
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system is a feasible design; however, aluminum-tubed systems 
could not provide a feasible design for the same output. 
Titanium tubed plants failed to produce a feasible design 
for a 30 MW (net) output power system. In both cases of 
infeasible design, the constraint which was consistently 
Violated was turbine internal efficiency, set at 90% for 
current state-of-the-art design. 

c. The energy conversion and efficiency of design 
of a turbine-generator has a major effect on the overall 
system performance as indicated in paragraph b above. 

d. The cost/KW output for titanium-tubed heat 
exchangers is one third the cost/KW output for aluminum-tubed 
heat exchangers using a 30-year life-cycle criterion, with a 
10% annual inflation rate and retubing at 10 and 20 year 
marks with AL-5052 tubing. 

e. Aluminum-tubed heat exchangers have larger tube 
bundle volumes, with volumetric differences between aluminum 
and titanium varying from 26.1 to 11.8% for evaporators and 
23.2 to 7.4% for condensers over the range of net power 
levels considered. In both cases volumetric differences 
Gumanish as the system's net electrical output increases to 
20 megawatts. 

f. COPES/CONMIN has provided optimum designs for 
each incremental output power level. By manipulating the 
specified design variables, subject to imposed constraints, 


COPES/CONMIN has created designs whose geometry and operating 
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parameters cannot be scaled on the basis of net power output 
(10 MW). Therefore, designs for component geometry at 
increasing power levels based upon such simplistic scaling 
@eeeeria Wil] not achieve an optimum design with respect to 


the cost objective function. 


B. RECOMMENDATIONS 
bh Evaluate additional objective functions including: 
a. Minimize heat exchanger volumes. 
b. Minimize parasitic power losses. 
c. Maximize thermodynamic efficiency. 
d. Maximize net electrical output. 

2. Perform a sensitivity analysis on power system design 
variables to evaluate their influence on component and system 
performance. This allows the designer to prioritize system 
components which can provide improvement in the objective 
function for a corresponding improvement in component design. 

3. Considerable uncertainties are associated with the 
expressions used to estimate component performances (two-phase 
pressure drops, film coefficients, etc.). The code should be 
tested to determine the sensitivity of system design to these 
uncertainties. 

4. Expand the code to include the use of enhanced heat 
transfer techniques and evaluate the influence of increased 


piping friction factors on pumping power requirements. 
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S. Evaluate proposed OTEC designs using proposed system 
parameter inputs, comparing both the basic analysis and the 
Seeimization output. 

6. Select other analytical expressions for heat transfer 
coefficients to validate the performance and output of the 
existing code. 

7. Evaluate the effect of a smaller thermal difference 
seen by the power system and its influence on a feasible 
design for a specific net electrical output. 

8. Evaluate the cost aspects of uSing variable-pitch 
pumps versus fixed-blade for a variable thermal gradient 
environment. 

9. Evaluate and verify the influence of incremental 
improvements (percent) in turbine internal/adiabatic efficiency 
with respect to gross and net electrical outputs and compare 


with the results reported in Ref. 16. 
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SEMemeelNEUT DATA FOR OTEC 


APPENDIX A 


EVAPORATCR - HORI ZUNTAL 


TUBE O.D. 
TUBE LENGTH 
SW TUBE VEL 
OPER PRESSURE 
TUBE MATERIAL - TI 
THERMAL CUOND(K) 


L.QUO( IM) 
+O.0U0(FT) 
6.Q0CO(FT/S) 


L30.Q000(LBF/IN2) 


TANTUM 


PeOOO(BTI/HR.FT.F) 


ANALYSIS 


2564 I0 0474) 
12.192 (ita 
1.829(M/S) 
0-830 (MPA) 


16.502(a/M.C) 


Meipes PROFILE = STAGGERED EQJI-LATERAL 


PITCH RATIO 

ENHANCEMCNT —- PLAl 
CONDENSER - HORIZONT 

TUBE O.D. 

TUBE LENGTH 

SW TUBE VEL 

OPER PRESSURE 

TUBE MATERIAL - TI 

THERMAL CONDO (K) 


1.30 
Ne UBE 
AL 
LeQUOCIN) 
56¢500(FT) 
6-000(FT/S) 


SI-OOO0(LBF/IWN2) 


TANIUM 


9-S300(BTU/HR.FT.F) 


252.490 (4M) 
17.221 (M) 
“LeS329{M/S) 
O-.514 (MPA) 


MOUBE PROFILE - STAGGERED EQUI-LATERAL 


PITCH RATIO 
ENHANCEMENT -— PLAT 
SALT WATER HOT PIPE 
Pare 1.0. 
Beerpe LENGTH 3 
SW PIPE VEL 
SW [INLET TEMP 
OW SALINITY 
SALT WATER COLD PIPE 
PIPE 1.D. 
wUPe ENG) H 
ow P[PE VEL 
SW INLET TEMP 
SW SALINITY 


pes0 
N TUBE 


L9.3900(FT) 
GO.-O00(FT) 
4.500(FT/S) 
BO.000(DE F) 
35. U/V00 


L13-600(F T) 


3000.000{(FT) 


De SVO(FT/S) 
40.000(DEu F) 
35. 0/000 


oC 


43 O3)(44) 
71.440 (id) 

Le372(M/S) 
266567(DEs C) 


50974 (i) 
314-400 074) 

beoroums 5) 

4.444(DEG C) 








AMMONIA CIRC PIPE 


PIPe I.0. aeJ00 Cr |) 

myree LENGTH L50.-000(FT) 
AMMONIA RE-FLUX PIPE 

mire 1.0. 2-000(FT) 

PIPE LENGTH DS0.0CO(FT) 


PUMP ANDO GEN-TURB PERFORMANCE 
EVAP Sw PUMP 
SeerclENCY MECH 3$5.00(PCT) 
COND Sw PUP 
Peer eENCY “eCH 85.00{PCT) 
AMMONIA CIRC PUMP 
SemecleNCY MECH 75. J0(PCT) 
Me 1URB EFFICIENCIES 
GEN MECHG&ELECT 96-650(PCT) 
TURB MECH $9. 80(PCT) 
POWER REQUIREMENTS 


0.610 (M) 
45.720(M) 


0-510(M) 


Loee +0 U4) 


ACC Feo re 1) 


AOTOR, Se soo GPCT) 


Gees Ore Gir Cella 


NET POWER QUTPUT 15.000( Mw) 


Ist 





ee ee 


oN x 5 


iPEE OTEG ANALYSIS OPTIMIZATION OUTPUT DATA 


EVAPCRATCR - HCRIZONT AL 


HT ABSORB 2905013760-0(3TU/HR) 587.2601 ( Mw) 

SW FLOW 3340 81230.-0(LBM/HR) 151535904.0(KG/HR) 
SW TEMP IN 80-Q0CO0(DEu F) 26.567(JEG C) 
Sw TEMP QUT 73e6724(DEG F) 23-i130(DEG C) 

NH3 FLOW 3788703-0(LBM/HR) L7L8519.I(KG/HR) 

OPER PRESSURE L30.0S7(LBF/IN2) 396293 7( KPA) 
EvVAP SAT TEMP ies o> OEG F) Z21-.435(0E% CC) 
PU eET FEAP T9-514(0EG F) 2leo7i (DEG NC) 
OUTLET QUALITY 92-00(PCT) 

NH3 PRESS DRuUP O-LS2(LBF/IN2 ) LeLISt( KPA) 

TUBE CHARACTERISTICS 
OUTTER DIA Q.952(1N) 24e1LB90 MM ) 
WALL THICK O.025(IN) 9.639(MM) 
LENGTH 42e132Z2(FT) b2aoo fi) 


MATERTAL - TITANIUM 
MIGeEsPROFILE = STAGGERED EQUI-LATERAL 


PITCH RATIO 1.+0 
ENHANCEMENT ~- PLAIN TUBE 
SW VELOCITY SsO2otrT/ Ss} begs (ey 
T WALL(SHELLSIDE) 71.456(DEG F) leg oowOEGnC) 
FILM TEMP (OSA IOUDES Fi) Z\eoGlit DEa i) 
DELTA T 3UILING Ue Jo2 DEG F) eZ) (VEGmwCy 
LeMeT.De de 7S57(DEG F) Sel9S( DES Ci 
EVAP EFFECTIVENESS 0.667 
NR OF TRANSFER UNITS 1.093 
OVL HT CUEF 612.97(B8TU/HR.FT2.F) 3480.58(W/42.C) 
H( WATER) Piero inher i 2s) S336.71(N/MN2.C) 
H( FOULING) Disa woewmereshneFIl2.-) 2L512.09(nN/12.C) 
H( METAL) 4409. 72(BTU/HReFT2.F) 25039.29( w/42.C) 


HC AMMONTA) 4093.4¢4¢(BTU/HR.FT2.F) 23215.00(W/A2.C) 


Fo2 





| 
| 





HT SURFACE Si aos DOLE T2 ) > selena ent 2M ) 


TUBE SHEET OIA laec Lor T ) 32294( 4) 
TOT NR GF TUBES 34449. 

Swear seESS JROP 4.06L(LBF/IN2) 23-00 5( KPA) 

MOISTURE SEPARATOR-INSIODE EVAP SHELL 

OPER PRESSURE L129.935(LBF/IN2) 395-872(KPA) 
QUTLET TEMP POeSS5S(0EG F) 21-296(0EG C) 
OUTLET QUALITY SJe 20(PCT) 
NH3 PRESS DROP O.4¢LlLO(LBF/IN2) 2e353( KPA) 


CONDENSER ~ HORIZONTAL 


enc ECT 1935472896.0(3TU/HR) 207.2210 Mn) 

SW FLOW 334871552.0(L34/HR) 15 1894363.0(KG/HR } 
Sw TEMP IN 4Q.000(0EG F) 42444(0EG C) 
Sw TEMP QUT 46.055(0EG F) Te308(DEG C) 

ie FLCW 3738708 -O(LBM/HR) 17138519.00KG/HR) 

OPER PRESSURE OF-LSL(LBFE/IN2) 607.7310 KPA) 
COND SAT TEMP 49.351L(DEG F) 9-0%5(DEG C) 
GUTEET TEMP 49.238(DEG F) 9-577(DEG C) 

NH3 PRESS OROP O-206(LBF/IN2 ) L.#23( KPA) 

TUBE CHARACTERISTICS 
QUTTER DIA 2e-972( IN) 2%-2533( 44) 
WALL THICK Q.025(IN) Je 542(0 M4) 
LENGTH Sart) Olt ile 17.500( 4) 


MATERIAL ~ TITANIUM 
TUBE PROFILE - STAGGERED EQUI-LATERAL 


PITCH RATIC 1.40 
ENHANCEMENT — PLAIN TUBE 
Sw VELOCITY 6-OLT(FT/S) 1.334( 4/58) 
T WALL(SHELLSIDE) 43.331L(0E5S F) S20 SCUEG CC} 
= ILM TEMP 438.785(DEG F) Peocolveo ©} 
weerA Tt COND UeZOT( DES F) O.504(DEG C) 
LM.T.0. 5e533(DEG F) Z3.L57(0ES C) 
COND EFFECTIVENESS 0.655 
NR OF TRANSFER JNITS 1.065 


5 








CVL HT COEF 

H( wATER ) 

H{ FOULING) 

H( METAL) 

H{ AMMONIA) 
HT SURFACE 
TUBE SHEET DIA 


TOT NR OF TUBES 


SW PRESS JROP 


446.-33(BTJ/HR.FT2Z.F) 
704.35(BTU/HR.FT2.F) 
Boece OUCBTIU/HR.«FI2.F } 
4393.63(8TU/HR.FI2.F) 
3033-55(B8TU/HR.FT2.F) 


762199.25{FT2) 


27-194(FT) 
Deli 2 
De S36(LBF/IN2Z) 


SALT WATER HOT PIPE 

PIPE 1.0. ZOsOTTET) 

PIPE LENGTH 300.000(FT) 

SW PIPE VEL 4eS9T(FT/S) 

SW FLOn 3340 81230-0(LBM/HR) 

SW INLET TEMP 80.000(DEG F) 

SW SALINITY 35. 0/000 

Sw PRESS DROP 0. 322(L3F/IN2) 
SALT WATER COLD PIPE 

PIPE I.D. 13.622(FT) 

PIPE LENGTH 3000.900(FT) 

SW PIPE VEL 5.334(FT/S) 

SW FLOW 334871552.0(L84/HR) 

SW INLET TEMP 40.Q000(D0ES F) 

SW SALINITY 35. 0/000 

SW PRESS DROP 0.508{LBF/IN2) 
AMMONIA CIRC PIPE 

PIPE 1.D. 2.00L(FT) 

PIPE LENGTH L50.0GO(FT) 

NH3 FLOW 3788708-0(LBY/HR) 


NH3 PRESS DROP 


15-033(LBF/IN2) 


AMMONIA RE-FLUX CIRC PIPE 


muPEs 3 2D. 

PIPE LENGTH 
NH3 FLOW 

NH3 PRESS DROP 


2.000(FT) 
59.0CO(FT) 


LL36612.0(LBM/HR) 


9-8 74(LBF/IN2) 


154 


2p ie eat 742 .C } 
BO seta (ts 42s C ) 
ZOO bet Ot N/ 12 .C) 
24947.99(A/M2-C) 
17367.903(W/42.C) 
79809.69( 42) 
3.239(%) 


3323611 KPA) 


Get J.) 


91.4491 '4) 
Le401L(M/5) 


151535904. 9( KG/HR) 


76000 1 CUE G 


2e21L7(KP4A) 


906/75( 4) 


914.-400{ 4) 
Le629(M4/S) 


151894 363840( KG/HR ) 


4et4e(DEG C) 


3eD01L( KPA) 


9-510(M) 
G5e729(M) 
L71383919.0(KG/HR ) 
1903.650(KPA) 


0-510(.4) 
152240 4) 
515555-3(KG/HR) 
68-079{( KPA) 








PUMP AND GEN-TURS PERFURMANCE 
EVAP SW PUMP 
HEAD. PRESS 2) x8) oie) Gaeta) 
war ACITY 653259.0(GAL/MIN) 
EEE CreENCY MECH 435.00(PCT) 
COND Sw PUMP 
HEAD PRESS 20) Tse at) 
CAPACITY GS211L9.2(GAL/MIN) 
Peeve leNCy MECH 35.00(PCT) 
AMMONIA CIRC PUMP 
HEAD PRESS PAze S521 rT) 
CAPACITY E2EI1. (GALS HAIN) 
mage elLENCY MECH %5.00(PCT) 
AMMONIA RE-FLUX PUMP 
HEAD PRESS 33-O015(FT) 
CAPACITY 5732 .64(GAL/MIN) 
BmeLCieNCy ™ECH 75.00(PCT) 
SEN-TURB EFFICIENCIES 


GEN MECHGELECT 96 .690(PCT) 
TURB MtCH VISGOCPC YT} 
TURB INTERNAL 39.83(PCT) 


TURB OUTLET QUALITY 96.77(PCT) 


POWER REQUIREMENTS 
TURB—-GEN GROSS 27663.313(HP) 
Srl ClENGyY LOSSES 


EVAP Sw PUMP 196%.851 (HP) 
COND SW PUMP #29 ~oS13(HP} 
NH3 CIRC PUMP 5941.714(HP) 


Neier e-ELUX PUAP 29.097(HP) 


NET POWER QUTPUT 


Coes Jt LIT /AIN?} 
MOTOR 98.090( PCT) 


Fe227(i4) 
2toge?i1le)( LET /MIN)} 
MOTOR 98.00( PCT) 


64.717(M) 
79505. 0( LIT/ MIN) 
MONO go. 00( PCT} 


L1L1.537() 
BS4IZ1<1(L1T/ MIN) 
MOTOR 98.99( PCT) 


20.633 (4wW) 
905594 Ma) 
1.496 (MW) 
3.143 (Mw) 
0.412 (AW) 
0.022(MwW) 


15.090 (MW) 


PERCENT PARASITIC POWER 24529 (PCT) 


THE MODYNAMEC CYCLE EFFICIENCY 


L355 


2e65(PCT) 








ost JF COAPCNENTS 


EVAPCRATOR 
CONDENSER 
GEN-TURBINE 
GENERATOR 

cVAP SW PUMP 
COND Sw PUMP 

NH3 CIRC PUMP 
NH3 RE-FLUX PUMP 


GPTIMUM COST 


Meat PER Net Kw OUTPUT 


8225672.90(0DCLLARS)} 
3567154.900 (DOLLARS) 
1578 775-00{D0LLARS) 
937205.06(DOLLARS) 
653332294(DOLLARS) 
§652298.06(DOLLAKS) 
136824.94{(D00LLARS) 

64443.34(DOLLARS) 


—=_——Vew es ee 


208 *+9216.00(D0LLARS) 
bo5 9. 75> ( 0UEL ARS) 


e310 





AEP ENDEX C 


SAMPLE COPES OPTIMIZATION AND SENSITIVITY ANALYSIS DATA 


$BLOCK_A_(TITLE CARD) | 
GECAN THERMAL ENERGY CONVERSION {OFEC) POWER SYSTEM 
$BLOCK 3 (PRUGRAM CONTROL PARAMETERS) 
2916,16 
2 16 L6 | 
$BLUCK C CINTEGVER OPT CONTROL PARAMETERS) 
5909075 
5 0 Q 5 
a OD {FLOATING PT OPT PROG PARAMETERS) 
0.9 
Q.V ae 
$BLOCK € (TOT NR DESIGN VAR,DESTGN O8J' [DENT AND SIGN) 
L6,27,- le O ‘ 
L6 a at -~1.0 - 
Mmeeen Ff (DESIGN VARIABLE SOUNDS,TNIT VALJES & SCALE FACTOR) 
1.0, 1.0420 | 
L.9 1.-Q+20 
1.0,1204+20 
1.0 L.0+20 
L.0,1.0+20 
hte, Ler cc 
1.9, 1.0+20 
le 1.0+20: 
85-02148.0 
85.0 143.U 
85.0%1148.0 
35.0 148.0 
ee 2 6 5 a B 
02D ZeD 
025,2.5 
0.5 2-5 
10.0,1.0+20 
10.0 ie Gio 
10.071.0+#+20 
LQ.V be O0t20 
2eJ,10.0 
20 10.0 
Zale 10.0 
220 10.0 
2-0,10.0 
2.0 LU.9 
2eJ,10.0 
2-0 10.0 
Letty 320 
1.+ 320 
1.4, 3.0 , 
1.4 aa) 
$BLOCK G (DESIGN VARIABLE IDENT) 
ly L,1.0 
a l 1.0 
292? r.0 
2 Z L.O 
3931.0 
3 3 L.90 
494e1.0 
; 4 4 1.0 
5932129 
5 5 L.O 


ad 








miro, Le O 


6 P) Leo 
fies 1 2 O 
i q 1.0 
miro, Le O 
3 3 UA: 
G99 1.0 
9 9 kes) 
10,10,1.0 
LO Lu 1.0 
im, i1,1.0 
Ll Ll 1.0 
L2,l21.0 
12 12 Lao 
Ladue L3, Le) 
L3 13 Law 
14,141.90 
14 14 LeU 
Powe o,y 1.0 7 
L 15 Va 
16,16,1.)0 ; 
lL 6 ALS, 
STRAINED PARAMETERS) 


y 

6 L 
Baeetk eink aoe CON 

> 

( 


— I. (CONSTRAINT [DENT AND BOUNDS) 
L? 
dew eae Bie 20-0 
l. OAS, Bie 0.0 
18,29 
13 20 
tee, 0. Or met Z20,0.0 
9-0 eit 2 2 0.V 
21523 
21 25 
0220.0 L~O0+2020-0 
Che 0.0 Ll. 0+290 0.0 
24 
2% 
30. 90-0,90.,7, 044) 
O. 0.0 SON 0.0 
22726 
a 
LOe s0202120+2070.0 
On. 0.90 1.9+29 sae 
BOLOCK P (SENSITIVITY OBJECTIVES) 
27,0 
atl Q P 
Meee ror 809+ a OG OID a i Do a 
9 LO Ll Ue 13 
bel 18 19 20 21 
25 26 DoT 
Pie 2223924925926 927 
fe) T 3 
14 L5 L6 
22 Zo 24 


$BLOCK Q (SENSITIVITY VARIABLE BOUNDS) 


9 


il 6 
LIS ee lO rlLSe.r1ldox 2009256 
io LO. 


15. lo. 20-6 
296 
é 6. 
18.,10. ldo: LOwe CU es 2D « 
lo. LO. 15. la. 20. 
3,6 


alESS 


2% 


Los 








3 fe) 
Pepe 2s Leslee x dep te — 
Zn QD le 
499 
4 6 
Wee oe De lerce 3 0c te _ 
Je 0.5 Ll. 
536 
5 6 
Moe s128.7129--130-r,131L- 2132. 
130. L23-e 129. 
6,5 
6 6 
89 .9871.2:8392259.,90er91L. 
89.6 87-6 83.6 
7,6 
7 6 
meg es (Dele ele 2501-5 
Le 05 0/5 
8,5 
; 3 5 
LevseDeeolSelevle25y1ledD 
le 2 Ve 
975 
9 6 
F Pees 9 3D0 9 400 9 4542 50~ 
4) 6 30. 35.6 
1095 
iD. 6 
D5. 94D 6 9206 9950¢7E007056 , 
55-6 45-6 52. 
11,6 
Ll 
meas Oe 7s leo ep Bes Se 
Te De oe 
12,6 
L2 6 
on 6 Zo oe >to ep Der Oe 
te ae a 
1346 
ies Dee Oe fie) Goto « 
7. De Oe 
1446 
oe 6 
pe este 8 Hoe Dee He 
ye ae Bee 
15,6 
| 15 5 
meOomwle tel eS ea leGelelyle 8 
1.6 1.4 toe) 
16.5 
| 16 6 
ee Ols 1.471.533 Le6y ie le I. 8 
1.6 1.4 Ee, 


END 





ae ein ya aie 
at OS = at O 
POO + 
PLPNOFSO 
ooo0oo 
=) =) =) 
aa aaa 
a2aee 


A 

S 

? 
NOMENCL ATURE 


YOO OU OVO 


Oe 6 oo Ss oo OOO OOO OC OOO OOOO O OOO oOo OoOoOoOoo }& 
CO CSP et’ OOO A AIUD OO POM SFO AIO OO NI ed CN et (OOO Nt Pe LA SOOO 
0 HO MO OM AO DD ev ed ed PP ed AD CFO OO TD et PD SLND MMOS OS 
FPP PPD et FS FOS a ed tA HUN AIO NID 0 tO et et et IO SICQULIAC 
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